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Electron energy distributions resulting from a continuous primary 

source at high energies in helium have been calculated. Contributions 

from both the primary electrons and subsequent secondaries are included. 

In the cases analyzed, the primary electrons are taken to be produced 

via the ionization of helium by heavy-charged particles from the 

nuclear reaction ^n^ + ,Li^ + ? He^ which represents an important 

means of using the neutron flux from a nuclear reactor to produce a 

plasma. The calculation covers an energy range from zero to about 1 KeV, 

the highest energy of the primary electrons associated with such a source. 

Electron-neutral collisions are assumed to be the dominate collision 

-5 

process. While this limits the fractional ionization to 10 , this region 

is of interest for application in areas such as radiation-pumped lasers 
and the buffer-gas region of the gaseous core nuclear reactor. 

Balance equations for a finite system are developed • for the ' slowing- 
down region, and the appropriate collision kernels are evaluated. These 
equations are solved numerically for the electron-flux distributions for 
primary source energies from 100 eV to 1 KeV, gas pressures. from 0.1 to 10 
Torr and also various plasma dimensions up to and including an infinite 
medium. Calculations of the W-value from these distributions for electrons 
in helium are found to be within 10% of experimental measurements. The 
results are also generally consistent with previous Monte Carlo calculation 
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A Boltzmann equation with appropriate source and sink terms is con- 
sidered for the low-energy range for cases both with and without a super- 
imposed electric field (0 <_ E/p <_ 10 V -cm * -Torr *). Without the field, 
major deviations from a Maxwellian distribution are only observed for low 

gas pressures ( < 0.4 torr), small systems (tube radius < 1 cm) and 

a, 13 -2 -1 

moderately high neutron flux values (> 5 x 10 cm -sec ). For the 

cases with an electric field, deviations from a Druyvestyen-type distribu- 

% 14 -2 

tion are observed when the neutron flux density is high (> 5 x 10 cm 

-2 .'v -1 -1 

-sec ) while the electric field is small (E/p< 4 V -cm Torr ). In 

general, the shifts of the distribution due to leakage and recombination 

depreciate the density in the low-energy region. This is attributed to 

the leakage of electrons during the slowing process. 

It is suggested that the methods developed here can be used to aid 

Monte Carlo calculations through the development of a combined analytic- 

Monte Carlo approach. The methods described can also be extended to 

applications involving different gases or gas mixtures, although the 

derivation of new collision kernels is required. 
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CHAPTER I 

DISCUSSION OF THE PROBLEM 

A. Introduction 

1 . The Problem and Objectives 

The present study is concerned with energy distribution of high 
energy electrons as they slow down and thermalize in a gaseous medium. The 
energy distribution in the entire energy range from source energies down is 
studied analytically. In particular, attention is given to a helium medium 
in which primary electrons are created by the passage of heavy-charged 
particles from nuclear reactions. As discussed later, a radiation-induced 
plasma of this sort is of interest in a variety of applications such as 
radiation pumped lasers and gaseous core nuclear reactors. 

A continuous source of high-energy electrons will establish a quasi- 
steady-state in the energy distribution of the entire electron population. 

The problem is complicated by the fact that this population not only includes 
the primary electrons but also secondaries created by ionization collisions 
during slowing down. This situation has been analysed not only for an 
isolated system but also for the important case where an external electric 
field is imposed upon the medium. 

2. Previous Work Related to the Problem 

Electron energy distributions in the absence of high-energy sources 
have been extensively analyzed^’ ^ ^ These studies all consider 

the situation where the distribution is sustained by an electric field, 
such as in the case of a gaseous discharge tube. On the other hand, the 
problem of electrons bom at a high energy has only been studied in the case 
of B-irradi'i: ^ns or for electrons created by gamma irradiation^ 
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However, in these cases the energy range considered has been restricted to 

quite high energies from MeV to the low keV region — since the interest 

has been mainly from the point of view of shielding calculations. 

Lacking complete electron energy distributions, previous workers 
interested in applications of radiation-induced plasmas have generally 
estimated electron ionization and excitation rates assuming a thermalized 
electron swarm with a characteristic temperature ^ ^ ^ ^ . While it is true 
that the slowing down time of electrons is frequently short so that a bulk 
of the electron population approaches a Maxwellian shape (no electric field), 
the tail of the distribution determines these reaction rates, and the assump- 
tion of a Maxwellian "tail" can lead to erroneous results. Stated another 
way, we note that electrons bom at high energies suffer inelastic collisions 
during the slowing process, and the resulting ionization and excitation rates 
are not adequately described by the Maxwellian tail associated with the ther- 
malized electrons. 

B. Physical Applications 

Plasmas resulting from electrons slowing down after creation at high 

energies are of interest in many recent applications. In a laser system 

which utilizes nuclear energy input, the primary electrons are created in the 

r oi 

lasing gas through the passage of heavy-charged particles through the gas LOJ . 

[■9I 

In the nuclear light-bulb reactor concept , the primary electrons enter as 

B-radiation or from the ionizations due to fission fragments passing through 

the fuel region and through a buffer gas region. In the electron-beam 

[ 10 ] 

fluorescence technique which is used to study basic parameters in gas 
dynamics, primary electrons are introduced into the system in the form of a 
beam from an electron accelerator or "gun". 
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The slowing down processes in these situations differ. In the laser, 
electron-neutral gas atom processes generally dominate. For the nuclear 
light-bulb reactor, the thermal temperature is high and the uranium is ther- 
mally ionized. This leads to a high degree of ionization and electron inter- 
actions with uranium ions play an important role. In the electron beam 
fluorescence technique, the spatial geometry and the space charge of the beam 
become important. However, in all of these systems the problem of high energy 
electron slowing down and the resultant total electron energy distribution is 
of common concern. 

The present study calculates the time and space independent electron 
energy spectrum for electrons slowing down in Helium. In particular, refer- 
ences are made to a system in which electrons are created by the passage and 
slowing of nuclear radiation as in a nuclear excited laser system. However, 
with some modifications, the present study can be applied to other situations. 
One situation of interest now at the University of Illinois is the electron- - 
energy distribution in the region close to the fuel boundary and the buffer 
gas region of a nuclear light -bulb reactor. 

C. The Model 

A variety of schemes have been suggested to use nuclear energy in 

r 8 ] 

achieving laser excitation 1 . One method employed by the University of 
Illinois group is to utilize the neutron flux from a reactor for the nuclear 
reaction n + £3^ ->• a + Li • Thus a laser tube is coated with Boron-10 

and inserted into a nuclear reactor. The resultant heavy-charged particles 
traverse the laser gas and deposit their energy which may enhance laser 
excitation. Over half of the energy deposited is in the form of ionized 
6-electrons. The resultant energy distribution of electrons studied here is 
therefore important to this concept. 
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: A continuous 6-electron source is assumed. Many of the actual 

experiments at Illinois have employed a reactor pulse. However, the full- 
width at half-maximum of the reactor pulse is of the order of ten milli- 
seconds whereas the collision time of electrons in helium at one torr is the 
order of micro-seconds or less. Thus the electron energy distribution at 
any time during the pulse can be viewed as a quasi steady-state distribution. 

While the present study primarily deals with the electron energy 
spectrum in helium, other gases or gas mixtures can be studied by the same 
method. Helium was chosen here because it is a main ingredient in many laser 
gas mixtures. Also it is representative of the noble gases such as neon that 
have been considered for use in the fuel and the buffer gas regions of the 
nuclear light -bulb reactor. Another reason for studying helium in this first 
analysis is that an abundance of cross section data available for it. This 
allows us to concentrate on the analytic technique, although even for helium 
great care is required in the development and selection of appropriate cross 
sections. 

Many of the physical parameters employed in the present study such 
as tube dimensions, the boron-coating thickness, and the neutron flux magni- 
tudes were selected to be consistent with radiation-induced plasma experiments 

r 12 1 r 131 

performed by J. C. Guyot et al. J,L , Other assumptions involved include 
uniform physical properties in the system and the approximation of the cylin- 
drical geometry by a slab geometry when energy deposition rates are computed 
{Figure G-l of Appendix G) . 

A Monte Carlo calculation has been performed using the same model ^ 


Both the present study and the Monte Carlo simulation serve to illustrate 
methods to obtain the solutions that can be applied to other situations of 
interest in the future. These methods compliment each other and depending on 
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the particular feature of a given problem, it may be that one of these 
methods is better from the point of view of efficiency and accuracy. For 
instance, in the high energy interval, the present analytical approach 
appears to be more efficient. This suggests that one future application of 
the present analytic method might be to combine it with the Monte Carlo code 
such that the high-energy region is treated analytically and this serves as 
a source to the low-energy Monte Carlo simulation. 
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CHAPTER II 

THE ELECTRON FLUX DISTRIBUTION AT ENERGIES ABOVE THE FIRST 
EXCITATION POTENTIAL OF HELIUM 


A. The Slowing-Down Flux Conservation Equations 

After creation via heavy charged-particle ionization or via direct 
introduction as a high energy electron beam in a gaseous medium, electrons 
lose their energy traversing the medium. The resulting distribution is 
described by the electron distribution function f (r,v,t) drdvdt, defined 
as the number of electrons in the phase space element, dr dv dt , at posi- 

— k _v 

tion r with velocity v and at time t. This function is governed by 
the Boltzmann Equation: 


It + * Ve^ 1 ) + I 


coll. (II-l) 


where F is the externally imposed force. The collision term on the right 
hand side describes the sources and sinks of new particles into the phase 
space due to collisions. These collisions include processes such as elastic 
scattering, excitation, ionization and recombination. 

To begin the problem at hand, we shall first consider a time and space 
independent case with no external force field imposed. In this limit, the 
energy distribution of electrons while slowing down is simply described by 
the collision term of the Boltzmann equation. The electrons are born at 
6-ray energies and suffer large energy losses through a series of inelastic 
-collisions. These include ionization collisions which produce secondary 
electrons that in turn slow down. Both primary and secondary electrons 
eventually thermalize, but we will first concentrate on the slowing down 
problem and therefore, we restrict the lowest energy to the first excita- 
tion potential of the gas. In this range, the electrons arc at a much 


"Page missing from available version" 
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included in the total flux calculation. Thus, instead of the approximation 
of equation (II-2), we consider the exact balance equation for collisions 
into and out of an energy interval dE at E due to the slowing down of 
electrons originating from, a delta function source at E . This is: 

Z tot (E)$(E)ae " s S(E-E 0 \U + £ z.CE'^ejSCE'jdt'de+SCEMlE (II ' 3) 

j J 

where, (E) = total collisional cross-section at energy E 

£• (E'->E)dE= differential cross-section for the scattering of 
electrons with energy E’ into energy dE at E through the j th type 
of collision. 

S(E)dE = source of secondary electrons introduced in the 
energy band dE at E. 

In practice it is most convenient to use an integral form of equation (II-3) . 
This eliminates the delta function and in effect smooths the cross-sections 
leading to less fluctuations in numerical integrations. 

From appendix (B) , the integral form of equation (II-3) takes the 


form: 


z\ SCEo^KtCE'jE) <Je‘ = S -f ( Ce„,E")K qe'Je^e" cn-4) 

\ "" c. -J i * P w 


where j denotes various excitation and ionization processes. 


Kj(E',E) = probability per track length that an electron suffers a 
jth type collision and slows down past energy E from E'. 

K s (E",E) •= probability per track length that a secondary electron is 
born with energy larger than E due to a collision by a primary 
electron with original energy E". 


The equivalence of this equation and equation (11-3) is demonstrated in 
appendix (B) . 
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Equation ( 1 1 - 4 ) states that the number of electrons slowing down past 

energy E must be equal to the introduction of electrons with energies larger 

r 16,17-1 

than E. Spencer and Fano and other authors J have worked with the same, 

type of balance equation for high energy 3-radiations (MeV slowing to keV) . 

They used the Holler 1 relativistic collision cross-sections and included 

Bremsstrahlung energy losses. In the present case, Bremsstrahlung can be 

neglected but energy losses from excitation collisions must be included. 

r i 9 201 

Also Vrien's Binary Collision Model cross -section L ’ , which is more 

appropriate for. the lower energy range, is employed. 

As shown below, elastic scatterings between neutral gas atoms and 
electrons give only a small contribution to the slowing down of electrons. 
[This effect, however, plays a role in the diffusion leakage of electrons 
out of a finite system, and will be considered later.] This elastic 
collision cross-section, as approximated by the Lent relation ^ ^ , is: 

_ 4 . UI-5) 

Q e i * 

where, Q. = total inelastic collision cross-section, 

’ u.n ’ 

Q , = total elastic collision cross-section, 
x el 

h = Plank's constant, 

E q = ionization potential of gas atom, 
a = first Bohr radius , 
and V = velocity of impinging electron. 

For an electron with energy above tens of eV in helium, Q. and Q . are of 
roughly the same order of magnitude. Since the mass ratio between electrons 
and the neutral atom is about 1/7300, the energy loss due to elastic colli- 
sions can be neglected. This fact has been further verified in Appendix (C) , 
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where the average fractiona] energy loss per track length due to elastic 
collisions with neutral atoms is shown to be small compared to that for 
inelastic collisions with atomic electrons. This comparison is done for 
the energy range above the excitation potential such that the Born approxi- 
mation is applicable for elastic collisions and the Moller cross-section 
is valid for electron-electron inelastic collisions. 

B. The Slowing-Down Kernels 

The binary encounter collision theory employed here follows the 

r 22 23] 

classical theory of atomic collisions developed by M. Gryzinski L ’ J with 

, . . .. . . . [ 24,25,26] • . „ 

later corrections and discussions by various authors . This type 

of cross-section was selected in preference to purely quantum mechanical 
cross-sections because their relatively simple form permits analytic 
integration. Another reason is their apparent good agreement with experi- 
mental measurements, especially in the lower energy range. 

In contrast to the Thomson cross-section, the Gryzinski cross-section 
accounts for the motion of atomic electrons in the target atom. The orig- 
inal Gryzinski model, however, is asymmetrical for the two interacting 
electrons. Thus, if both the incident electron and the atomic electron are 

at a distance r from the nucleus, one would have zero and the other would 

2 [27] f 291 

have -e /r potential energy. To correct this, Burgess 1 3 and Vriens 1 

introduced a binary encounter model with symmetrical treatment of the two 
interacting electrons. Thus, the incident electron gains a kinetic energy W 
and at the same time loses the same amount of potential energy as it inter- 
acts with the atomic electron. The atomic electron is assumed to be bounded 
with the same potential energy, IV; thus the two electrons are symmetrical 
and indistinguishable upon interchange. 
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In a later refinement, Vriens^^ developed a quantum mechanical formulae 
for transfer of momentum and energy between a pair of free' colliding elec- 
trons. He then applied these formulae to the electron-atom binary encounter 
theory to obtain differential and total cross-sections. The symmetrical 
collision model was also incorporated and the influence of the nucleus was 
included. The resultant differential cross-section contains the direct, 
exchange, and interference contributions, and it is given by Appendix (D) . 
The differential cross-section for a loss of energy AE is given as: 


_ 1 + ^ 3 - 

C AE 4 °AE + 4 C ££ 


(H-6) 


and 


± .... fte 

O'.. d4E = — - 
At E, 


+ 4 


4 


Ue a £ 3 (e 3 -E 2 -ae) z 


-V i. 


2 §' , . 

i 4E(f 3 'E z -Ae) - E * 


3 


4 <UE = *g* 

E e 3 4 E z 




+ 4 ^ 

+ 3 


4E 5 CE 3 -E i -AE) Z 


CII-7) 


4_ 

3 


Z$' 






C £ -E 2 -4E) 3 lEs-Hi-AElAE: JV E 2 


S ')***! 


( 1 1 - 8) 


where 


E, = 


E„ = 


energy of incident electron 

average kinetic energy of atomic electron 


E, = 


E 1 + E 2 + U 


U = ionization potential energy 


and 




<5>' = 


E 3 - AE 


1/2 


cos 


E,-E_ 
o 2 


In 


E„-E - E 
o 2 


R = Ryberg's constant. 
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The total ionization cross-section derived from this differential 
cross section has been compared with experimental results, and good agree- 
ment is observed at low to moderate energies. However, in the present 
work, we are more interested in a direct check of the differential cross 
sections [Equations (II-7) and ( 1 1 - 8) ] since they lead to the slowing down 
kernels. Since direct measurements of the cross sections are not available, 
it was decided to compare calculated and measured energy distributions for 
secondary electrons produced by a primary beam of electrons of energy E on 
a suitable thin target. Figures ( 1 to 5) compare the calculated 
probabilities using binary models with experimental measurements by C. B. 
Opal et.al.^ The experiments were performed at a pressure of 2 x 10 ^ 

torr or less and multiple scatterings would be expected to be negligible. 
Quite good agreement is observed, lending confidence to our use of these 
cross-sections . 

Appendix D gives the calculation of slowing -past kernels by 
integrating the differential energy loss cross-sections over the appropriate 
limits. The probability per unit track length that an electron will emerge 
with energy less than E from energy E' due to an ionization collision, 

Kj (E ' , E) , is given by: 

for (E ' + U)/2 < E'-E; Kj(E',E)=0 (n _ 9) 

for (E* + U)/2 > E'-E > U; 


KjCE'.E) 


_ NeTie 1 


_j — b r - 1 

tE'-£ 3 1 L (E'-E,)' 


( E 3-E 


y] 


V 

(e’+U) 


In 


C-S3B} 


( 1 1 - 1 o ) 


and for (E'-E) < U: 



(cm 2 /eV) Cross-Section (cm 2 /eV) 
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Ejected Electron Energy (eV) Ejected Electron Energy (eV) 


Figure 1, Figure 2. 



Ejected Electron Energy (sV) Ejected Electron Energy (eV) 

Figure 3. Figure 4. 


Figures 1, to 4 . Comparison of Calculated and Pleasured Secondary- 

Electron Distributions for Primary Energies of 
100, 300, 500 and 1000 eV 
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E 3 


CE'+U) 


{l 3 " E' + 3 [tf E ,2 1 

4 •■ (f)) , 


(11-11) 


where we have used the same notations as in equations (II-7) and (II-8) with 

the addition that N is the total number of valance atomic electrons per cm" 3 . 

e r 

In the case of helium, N g = 2N where N is the density of helium atoms, and 
one does not have to sum over the electrons in different atomic shells as 
in the case of more complex atoms discussed in Appendix D. Also is 
defined in Appendix D. 

The probability per unit track length that an electron will emerge with 
energy less than E from energy E' due to an excitation collision is 
K ex (E',E). Within the one excitation level model discussed in Appendix D, 
the incident electron loses an amount of energy x which is larger than the 
first excitation energy but less than the ionization energy of the atom. 

The atomic electron does not gain sufficient energy to be ionized, an exci- 
tation is considered to have occurred, and the excitation slowing-past kernel 


is given by: 


for E'-E > U ; K ex (E',E) = 0 


for U > E'-E >_ U x ; 


(11-12) 


( _J 1 ' J ) 2 E r \ 

lE'-E U f 3 Z [ (E ,_ E yz 

_ J t I V- , n r U(etQ) - j 1 

U 7 (ew) z j ^ l+u ) Le' (E‘-E) J / 


( 11 - 13 ) 


and for >_ E' -E ; 

^ Ur r) -lir 


i i r 

u’+u-u.y- j t + u L 


U(.e*+ o- U,) 

£' U, 


( 11 - 14 ) 
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. Ml 

where is the first excitation level energy and $ is defined in 
Appendix D. 

The secondary production kernel is defined as the probability per unit 
track length that an electron with energy E' will produce, through an 
ionization collision, two electrons with energies both larger than E. 

This is given by Appendix D as : 


for 

for 


£' £ 2.E f U 
E 1 ^ ZE +-U 

k 5 (e!E) = 



Ce'.E) 


o 


Cl 1-15) 



When equation (I 1-4) is rewritten with the three slowing down past 
kernels considered above, we obtain: 


6 K I u , ,6)$tE')a E ' + \ K„ce:E)5<e')<U' = yJ 


■ t-*4 


E + U 




ZE*4) 

(11-17) 


where ^ = £ + VJ for 2.E U ^ E 0 and A ■= E*— E if 

2E + U^E 0 

The form of K can stand some improvement. When experimental 
e*. 

excitation cross-sections are available ^CtjeJcan be expressed exactly 
according to the experimental results. We shall let K^Ce^E^ be the 
probability that an electron will slow down past E from E' per unit track 
length due to an excitation collision which leaves the atomic electron 
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in the excited state. Then the relation between Ke* ^ £ ^ and 
tvl e O-i ( E 1 ^ , the total excitation cross-section at energy E' for the 
jth excitation level, is: 


Kt E 1 , E) = O 
e*. ' 


for E' « U 




KI«(e',E)= ^eO-'CE') W Ce’-E) $ U ^ ,E'^U< 
K^Ce'.e'l-o W (e'-e)^O^ 

_ -tk 


(11-18) 

(11-19) 

( 11 - 20 ) 


where VJ; — ^ 


level excitation energy. 

We are able to define K ft xC e ’, E) in this simple manner because 
a definite amount of energy is lost in an excitation collision. No such 
relation exists, however, for ionization collisions so the experimental 
total ionization cross-section cannot be utilized in this manner. With the 
notation of Ce’, E) , equation (11-17) takes the form: 

j KjCE', E)§CE‘)cJe' + \ K e>c ( ^Ce'^cIe' = 


Sp + \ Vc 3 Ce i ,e)^Ce‘He' 


(11-21) 

For the case of helium, we shall include the excitation levels with the 
largest cross-sections in our calculations. Although the experimental 
observations are incomplete, especially for n = 2, we shall approximate 
them with analytic expressions. We shall demonstrate that this is superior 
to the one-level Gryzinski-type excitation model. 


C. W 7 -Value and Leakage 

The production of secondary electrons and the excitation of atoms occur 
at an energy above the first excitation potential in neutral ground state 
noble gases. Thus the energy flux spectrum obtained through equation (1 1-4) 
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or equation (11-21) can be used to calculate the primary input energy 
necessary to produce one ion-pair. This is best done through %Ce = U.) / 

the slowing down past density (total # slowing down past energy E/time 
volume), which is ' as shown in Appendix E: 

% CE = 0,-) = (' $CE")kmCE", U.HW' (n _ 23) 

«J O i 

where is defined as the number of electrons that appear with 

energy below H per unit track length at energy E". With this, the energy 
necessary to produce an ion-pair, W, can be calculated as : 


VI — s 'i >e °/ /£ |)-Ce 


(11-24) 


The W-value found in this manner provides a convenient check with experi- 
mental data since a number of such measurements have been reported. 

Another aspect of the present problem is that leakage must be considered 
if the system is finite, i.e. especially when the system is small and/or pressure 
low., The inclusion of leakage losses with a spatially dependent source is a 
complex problem. We shall assume cylindrical geometry with uniform electron 

production throughout the gas volume. Then if leakage is approximated by the 

- , [12,131 

first fundamental mode as postulated by J . C. Guyot 1 in his atomic 

metastable density calculations, the spatial flux $Ct, is given by: 




L. / (11-25) ■ 

where R and L are the radius and length of tube and the origin of the 
co-ordinate system is at the center of the tube, and A is an normalization 
which depends on the source strength. 

The diffusion approximation of the particle transport equation 
(Appendix F) gives a diffusion coefficient in terms of collision cross-sections. 
For the high energies involved in slowing down, this diffusion coefficient 
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should be used instead of the ambipolar diffusion coefficient . Then the leakage 
term is approximated by: 

~ D(eW 2 5(Ec,E) — ^,CE a . E.') (11-26) 

where TV is a characteristic length of the system. For a cylindrical tube, 
within the fundamental mode approximation, is given by: 

= + (x) (11-27) 

With this leakage term included, equation (11-21) becomes: 

; 4 A __ 

e)$(e.,eI) +Z 

5 


(11-28) 

where the first term on the L.H.S. accounts for rate of leakage out of the 
system by electrons with energy larger than E. 



iE + U 





D. Computation and Results 
1. Infinite Medium 

The case when leakage can be neglected is considered first. 
Following Fano and Spencer equation (11-17) was solved numerically. The 
integrals are expanded by summations and equation (11-17) becomes: 


k‘ 

k:CEi En^CEt E^WiAE,- 




P=2*+I 


\ l (11-29) 


where are weights of numerical integration, T = gives the 

number of intervals in the ionization potential, and the index j denotes the 
jth type of slowing down process. The quantity k is equal to 2.n-'-T 
if zn + I^k ; or equal to k if a.n+1 > k ; f : being the maximum 6-ray 
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energy. For the last tern on the right hand side of equation (11-29) 

vanishes. Equation (11-29) is solved step-wise on a computer from n = k 
down. For a general W , is given by: 




r k h 

_ {.sm+I \ ; = J -J 


C^ n ,E n ))Vl A AE n 


(11-30) 


The flux density distribution of equation (11-30) has been solved using 
several different models of cross-sections in helium. These cross sections 
are illustrated in the following figures and table. Figure (6) compares the 
total ionisation cross-section obtained through the integration of Vrien's 
differential cross-section (binary collision theory) with experimental 
observation. A semi-empirical "modified" Vriens cross-section is also shown 
which uses multiplicative fitting parameters to match the experimental 
results. 


For excitation losses, since we can utilize the total experimental 
cross-sections, empirical formulas are used to fit the experimental data. 
Table ( 1 1 - 1 ) and Figures (7-10) illustrate some of the excitation cross- 
sectional values we used. We have only included the six largest excitation 
cross-sections as the others will give comparatively little contribution. 

When experimental measurements are not available, we have resorted to a 

Bethe-Born type approximation ^^>3l] _ 

Figure ( 1 1 - 1 1 ) compares the electron flux density spectrum for a unit 
primary electron source at 500 eV obtained by using a one excitation level 
model with that obtained using the excitation cross-sections .of Table ( 1 1 - 1 ) . 
The general shape of distribution is similar to those obtained by Miller^^l 
for 3-radiation in water, and it is also in general agreement with Monte- 


"Page missing from available version" 
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Table II-l. Fitted Excitation Cross-Section for Helium 


Excited State 

Cross-section (A 0 )^ 

Energy Range (eV) 

Reference for 
Data Used 
in Fit 

2 3 S 

0 . 0306 (E-19 . 8) A nr . 0 , 

0 , 2 s dFtOP* °' 0086 

Above 19.7 

[34] ,[35] 

2 1 S 

3.5 .. .... 0.0317 , 

E/R ! °-° 455 - F/R ) 

Above 20 . 6 

[30], [31] 

2 3 P 

0.0114(0-20.8) „ 0032 

0.25 + (E-20 . 8) 2 

Above 21.2 

[36] 

2 X P 

- 0.085+0.177 Jin | + 

0.0376 AU r 

— Above 21.5 

[30], [31] 

3 3 S 

[9.6E - 224] x 10' 4 

22.6-35 

[33] 


[-4.66E + 273] x 10~ 4 

35 - .65 



[2 x 10‘ 4 ] 

Above 65 


... 3 X S 

[3.4E - 78] x 10" 4 

22.8 - 37 

[33] 


[-0.45E + 64.8] x 10~ 4 

37 - 70 



[-0.112E + 40.8] x 10' 4 

70 - 207 



[-0.0208E + 22.17] x 10" 4 

207 - 1000 


3 J P 

[6E - 100] x 10" 4 

23.3 - 65 

[33] 


[1.1E + 200] x 10" 4 

65 - 100 



[-0.433E + 373.3] x 10' 4 

Above 100 



cr(cm 2 ) x 10 18 a- (cm 2 ) x 10 10 
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f 15] 

Carlo calculation by B. Wang L . i'he rise of the flux density after about 

150 eV is entirely due to the creation and slowing down of the secondary 

electrons, W-values have been calculated from these distributions using 

the slowing-past kernel developed in Appendix E. The W-value 

obtained from the distribution from the experimental excitation' cross-section 

is 43.5 eV per ion pair and the one obtained from the one excitation level 

model is 49.1 eV'/ion pair. This compares with the experimental value of Jesse 

r 37] 

and Sadanskis of 42.3 eV/ion pair L . This demonstrates that the one 
excitation. level model , which is admittedly crude, tends to overestimate 
the amount of excitation. 

Figures 12 and 13 show the flux density distributions for a 1 Kev and 
a 250 eV primary electron source. As indicated in table (II-2), the W-value 
from 250 to 1000 eV varies less than 10%. This verifies experimental indi- 
cations that the W-values are essentially energy independent over this 

r 38l 

range 1 . Furthermore, since the ratio of excitation cross-section to 
ionization cross-section is independent of pressure, the W-value will be 
pressure independent in this model. 

Figure 14 compares the results obtained using the binary collision theory 
cross-section (Vriens) with those for the modified Vriens or semi-empirical 
fitted cross-section for a primary source energy of 1 keV. Figure 15 
compares the above two cases together with a one excitation level model case 
for a source energy of 500 eV. 

The curves for the two models shown in Figures 14 and 15 are reasonably 
consistent, although some discrepancies are noted. Some further feeling for 
the agreement of the two can be obtained by calculating W-values from these 
distributions. The W-values are listed in Table II-2. Both results are 
close to the experimental value and are essentially independent of energy. 





ELECTRON FLUX DENSITY (cm sec' eV ) 
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ENERGY (eV) 

Figure 13. Electron Flux Distribution for a 250-eV Primary Source 
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Exp. Excit X-sect. 

Gryzinski 1 Level Excit. 
Modified Vriens with Exp. Excit. 
E. = 500 eV 


ENERGY (eV) 


Figure 15, Comparison, of Electron Flux Distributions Obtained 

Using Multi-level Experimental Excitation Cross Sections 
with 1) Vriens and 2) Fitted Vriens Cross Sections and 
by a 1-level Excitation Model 
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Table 

I I -2 . W-Values for Electrons in 

Helium 3 

Primary energy 
(eV) 

Calculated W-value, eV/ion pair 
Vrien's cross section 

Semi -empirical cross- 
section 

1,000 

41 (49. l) b 

44 

500 

43.5 , 

45 

250 

44.5 

45 


a The average W-value over this energy range as measured by 
Jesse and Sadankis is 42.3 eV/pr [37] 

bW-value calculated with the one excitation level model. 


Table II-3. W-Values for Various Noble Gases C 

W-value, eV/ion pair 

Gas Ionization Energy 1st Excitation ?Energy Calculated Experiment L 



He 

24.47 

19.7 

49.1 

42.3 

Ne 

21 

15.7 

42.9 

36.6 

A 

15 

11.2 

35.6 

26.4 

Kr 

13.99 

10.0 

31.2 

24.1 

Xe 

12.13 

9.0 

25.6 

22 


c Primary electron energy of 1 keV 
^One excitation level model 
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Thus it can be concluded that either cross section model can be used with 
reasonable confidence. In contrast, as observed earlier, W- values for the 
simple one-level model are considerably higher than the experimental values. 

Figure 16 shows the electron flux density distribution for three initial 
energies: 1 keV, 0.5 keV, and 0.25 keV. These curves were obtained by using 

the six excitation cross-sections and semi-empirical- form of the binary 
encounter collision cross-section. As expected the shapes of the curves 
are similar but shifted in proportion to the source energies. 

As stressed earlier, we view the W-value as a means of checking the 
validity of the flux distribution calculation. As a matter of interest, 
this has been extended to noble gases other than helium. Since excitation 
cross-sections are incomplete, the one-excitation level model and the 
Gryzinski energy loss cross -section^" have been used in all cases. Table 
(II-3) compares the results with experimental measurements. Since, as 
demonstrated earlier for helium, the one- --level model is somewhat inaccurate, 
the results are only expected to display trend%. Indeed the calculated 
W-values are consistently high, but the rough agreement for this variety 
of gases gives confidence to the present treatment. 

Figure 17 compares the present calculation for a 1000 eV primary elec- 
tron distribution with some normalized results from a Monte Carlo simula- 
tion 1-15], p,| onte carlo technique may be regarded as a mathematical 

experiment and the result is presented as points in Figure 17 ^ The -general 
agreement in shape is encouraging, and the discrepancies are thought to be 
due to the rather simplified cross-sections employed in the Monte Carlo 
code' 15 *. ' 


E 0 = 250 eV 



E 0 = 1000 eV 


With Modified Vriens 
Cross-sections 


E 0 = 500 eV 


ENERGY (eV) 


1000 


Figure 16. Electron Energy Distributions Calculated with 

Experimental Excitation and Fitted Vriens Cross 
Sections for 250, 500 and 1-KeV Primaries 
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KX) IOOC 4000 

Energy (eV) 

Figure 17. Comparison of Electron Flux Distributions from a Monte Carlo 


Calculation and the Present Calculation for 1-KeV Primary 
Source 
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2. Inclusion of Leakage 

For finite media where leakage effects must be included, equation 
(11-28) is employed rather than equation (11-17). The numerical integration 
is similar however. In this case the total elastic cross-section is re- 
quired and it is estimated using the first Born Approximation. As derived 

[39 7 

from Appendix C starting from the result presented by Mott and Massey 1 , 


we obtain : 


3K 4 \’-CX , '^4k i ) 3 ~~ 


where k = rn„vJeA; 7 X = Z'z/cX.a <Xo = Bohr Radius. Values of this 

cross-section and some experimental measurements are compared in Figure 18 
for energies above 10 eV. Then for a cylindrical system of radius 1.27 cm 
and length 122 cm, \//t is estimated to be 3.6 cm -2 For this dimen- 

sion Figure 19 compares the flux 'distributions obtained with leakage to that 
for an infinite medium with a 500 eV primary electron and 2 torr pressure. 
Figure 20 shows the case for a pressure of 4 torr. Comparing Figures 19 
and 20 one sees that with increased pressure a finite system approaches one 
where leakage can be neglected. This is because a t higher pressures electrons 
suffer more collisions before they can leak out of the system. 

Figure 21 shows results for a primary electron energy of 800 eV case 
and a pressure of 2 torr. Comparison of Figures 19 and 21 shows that the 
diffusion loss is more significant for a higher energy primary source. The 
physical reason is that the higher energy electrons and their secondaries 
have a better chance of leaking out of the system. Thus for higher energy 
primaries one must be concerned about the leakage up to a higher pressure. 

Changing ‘/A. 1 is equivalent to changing the size of the system. 

Figure 22 shows results for three values of for the case of E q = 800 eV 


Electron Cross Section (X 2 ) 


1.5 



□ La Bahn and Callaway [40 ] 
Also L. Vriens el ol [40] 

O D.E. Golden etal [41 ] 

A C. E. Normand §2 ] 

— Born Approx. 


Figure 18, 
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Comparison of Calculated Total Electron Elastic Scattering Cross Section in Helium with. 


Measurements 
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• Figure 20. Comparison of Flux Distributions With and Without Leakage 

at p= 4 Torr, E = 500 eV 
r o 




Electron 
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Figure 21. 
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Comparison of Flux Distributions Vith and Without Leakage 
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10 100 1000 

Energy (eV) 

Figure 22. Comparison of Flux Distributions With E q = 800 eV, 
p = 4 Torr for Three Cases of System Dimensions; 
l^ 2 = 0, 3.6 and 7.2 cm -2 
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at 4 torr. A decrease in physical dimensions increases the leakage and 
shifts the flux to lower energies. The cases of ‘/a x = 0 , 3.6 and 7.2 cm ^ 
correspond to the infinite medium, a tube radius of 1.27 C m and a tube 
radius of 0.88 cm, respectively. 

These results are to be compared with that of Figure 23 for the case 
of 250 eV primary electron energy. (The gas pressure in Figure 23 is 2 torr 
because curves for higher pressures would overlap the infinite medium case.) 
A comparison of Figures 22 and 23 shows that the change of size of system 
has less effect on the energy spectrum in the case of a lower primary elec- 
tron energy. Another observation is that the leakage term has a cumulative 
effect and the percentage deviation between the infinite and finite media 
cases becomes more significant at lower energies. For instance, for the 

250 eV primary electron case, at 125 eV the difference between Zl't — O 

\ 

and Y& - 3.6 cm” is about 22% while at 50 eV the difference is about 30%. 
This is because the population is progressively reduced as electrons leak 
out of the system during the slowing process. 



Energy (eV) 

Figure 23. Comparison of Flu.x Distributions With E q = 250 eV 

p= 2 Torr for Three Cases of Sy ra Dimensions; 

2 - n 
1/a = 0, 3.6 and 7.2 cm 
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CHAPTER III 

THE ELECTRON ENERGY DISTRIBUTION BELOW THE FIRST 
-- ' EXCITATION POTENTIAL ENERGY OF HELIUM 

A. A Model for Thermal iz at ion in Helium 

In a weakly ionized plasma, one can expect to neglect electron-electron 
and electron-ion interactions during thermalization. The validity of this 
assumption will be examined later in this chapter for the specific physical 
situations of interest. In conventional gaseous discharges, electron- 
electron interactions are thought to play an important role in bringing the 
energy distribution closer to a Maxwellian. However, in the present case, 
it will be shown that the absence of an electric field allows the distribution 
to collapse into a nearly Maxwellian form. 

Due to the absepce of ionization and excitation collisions, below the 
first excitation energy, the electrons in an energy interval are populated 
mainly by elastic collisions and depopulated by recombination with ions, 
leakage and scattering out of that energy interval. The steady state equation 
taking into account of these processes is: 

•+ 2: T (E^§C£\« ^ r s (E , ^E)$(£'yat'-I s (e)$CE) (HI-1) 

o . 

where D(E) = diffusion coefficient 

<&(E) = electron flux distribution 
I r (E) = macroscopic electron-ion recomfoination in cross-section 
Z s (E) = total macroscopic scattering cznoss -section 
E s (E' E) = macroscopic scattering cross-section of going into 
interval E from E' . 

S(E) = direct electron source from high energy inelastic 
scattering. 
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In equation (III-l), the diffusion approximation notation has been used 
(see Appendix D) . Because of the small amount of energy less per elastic 
collision it is difficult to solve equation (III-l) by numerical integration 
techniques as was done in the high energy region. Instead, one can make use 
of the fact that energy loss per collision is small and expand $(E) in 
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moments of (E'-E). This procedure is similar to the "heavy gas model" 1 
in neutron thermalization. It should be noted that the assumption of small 
energy loss per collision breaks down when electron-electron or electron-ion 
collisions become significant. This model is, therefore, valid only in a weak 
ly ionized plasma. 

In the absence of leakage and recombination, the principle of detail 
(44] 

balance 1 - J shows that $(E) = M(E), where M(E) is a Maxwellian distribution. 
Thus for the case where there is recombination and leakage, it is convenient 

i 

to define $(E) as M(E)^(E) and solve for the correction factor i/>(E) • With 
this substitution equation (III-l) becomes: 

J (III-2) 

' c°° 

Normalization is such that \ where n is the total primary 

o 

and secondary electron density. Electrons are assumed to be born at higher 
energies so electrons enter the low energy region by slowing and an explicit 
source is not included in (III-2).if a low energy source term were included, 
equation (III-2) would be inhomogeneous. However, the solution for the 
inhomogeneous case can be obtained from the solutions of the present homo- 
geneous equation by the method of variation of paraiseters as shown in 
Appendix G. 

Applying the principle of detail balance to the Maxwellian factor of 
the distribution, we obtain: 
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, -E'/kT „ , 
Ee x 3 Ce‘ 


•E'i = E 


e '*X 


(E-*E‘) 


(III-3) 


The integral on the R.H.S. of equation (III-2) becomes: 

r 6 ' 

j I s C£‘-^p-')H j Ce")tAC&'')cXE' = < Ae.' (n'i-4) 

° • 

and substitution back into equation (III-2) gives: 


+r r Ct^CE^ = ^H s CE-'E l ms'}c)lE' -l^E^CE) f II 1-5) 


Expanding ip(E') around E’ = E in the integral of equation (III-5) by a 


Taylor series, one obtains 




(II 1 - 6 ) 


With this expansion, the integral on the R.H.S. of equation (HI-5) becomes: 


Jz,(E-.e’)ifie‘) 4 E* \z s ^i) + L ae«Z,Ce> jp) 


4*#<a to” 


n l 


(HI- 7) 


where, 4 E" = ^ 

The assumption next used is that the energy loss due to elastic scattering 
between electrons and neutral atoms is small. This is reasonable despite a 
large scattering cross section due to the large mass difference. The 

energy E’ is therefore close to energy E. 

m . 

0 

For jj = — << 1. where M is the mass of gas atom and m is the mass of 
M e 

r 43i 

electron and assuming isotropic scattering, the moments of AE are given by 1 J 



AE 1 = Ifx UT - 

E'l + & CfX 1 ) 


and 

AE 1 = 4 jx £ 

+ &(J*) 

- 

4E n < ere 


(II 1-8) 


^ (E-£) n <k' 


45 


where T is the gas temperature. 

In the case of elastic collisions between electrons and neutral atoms, the 

2 

terms with order larger than u are not retained, and equation (I I I -7) 
becomes : , 

= I s (6')4'fE'l +X J (E)[2UO.kT- E\)E <f( e ) 

+ U S (E (in-9) 

Also equation (III-5) becomes: 


-3>(e^ v SJ'Ce') = -1 t <;e)H>Ce) + } ZM) lUWT-E) ^(e) -+ CeY] 




(III-10) 


Changing the notation back to 4(E), where 4(E) = i|j(E)M(E), Appendix H 
shows that equation (III-10) can be written as: 


-KEl^SCEl-tl^Eme) = JZ^CElU ?Ce)ve 4 §CE)^-EfeT^-5^1] 


(III-1I) 


To proceed in solving equation (III-ll), we must assign specific forms 

to the cross-sections and -the diffusion coefficient. In the case of helium, 

the elastic scattering cross-section is reasonably constant in this low 
(451 

energy range . If we assume a constant recombination coefficient over the 
thermal range, the recombination collision cross-section i- s .. = CV v elocity) > i> e 
1/u. Although a weakly ionized gas is assumed throughout this development, 
for the low energy portion, it can be assumed that diffusion leakage is 
governed by the ambipolar diffusion coefficient. - 

One must be cautious that some of these assumptions may not be entirely 
valid in certain pressure ranges or if the model is applied to gases other 
than helium. Also, as examined in detail in later sections, the degree of 
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ionization must be below a certain critical value for equation (III-ll) to 
be valid. However, in general, when electron neutral atom scattering 
dominates, and when the moments of averaged energy loss can be cut-off at 
n = 2, equation (III-ll) is accurate. This equation can be applied to gases 
other than helium by use of appropiate cross-sections and coefficients. If 
an analytic solution cannot be found, equation (III-ll) can still be solved 
numerically. This is not necessary for helium, however, since with some 
simplifying assumptions, an analytic solution is possible. While the assump- 
tion may lead to some inaccuracy, this technique is preferred since consider- 
able physical insight is obtained. 


B. Series Solution of the Thermal ization Equation 

With the assumptions that the total elastic scattering cross-section 
is a constant and that the recombination cross-section is an inverse function 
of the velocity, as shown in Appendix H, equation (III-ll) becomes: 


- C I - z-tfl + r -4'd-V}' A l'£ < V = ° 

where 

T 


4= Mjg) 

3 2-s 



Here D is the ambioolar diffusion coefficient, and A is the characteristic 
a * 

first fundamental mode diffusion length defined in Appendix F. Equation 
(III-12) is solved by expanding 4>(y) in polynomials of y, i.e.; 


= l 

sofl - 2. N y 1 

iru 


where , 


(III-13) 
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The Maxwellian distribution 's used as a multiplier because for g = A = 0, 
i.e. no leakage and absorption, the solution of equation (III-12) is a 
Maxwellian as expected. Appendix (H) gives the . coefficients of the series 
S(y) to be: •/ 


and 



= A A «/3 


(III-14) 


Appendix (H) gives the first four coefficients' to illustrate the form of A^. 
It is noted that terms with higher orders of A and g appears later in the 
series. The assumption of weak absorption and leakage is not a necessary- 
one for the series to converge. However, for the physical cases considered, 
g and A are muchuless than one. In these cases, one can expand S(y) in A 
and g and retain the terms of first few orders in g and A. The series S(y) 
for keeping terms up to second order in g and A is : 




(III-15) 


The series ^(y), JU 2 (y) * ^ i Cy) > *[2^ and X(y) rapidly converge for y < 5 
and they have been derived in Appendix (H) . The final form of ^(y) is given 
as : 


= M* e l -h ^ 


,in+ I 


Un+o r U-* Vi) 



1 

(zb.+\)r(n-<- s '/i) 






(III-16) 
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where, b = '/% , b, n = ( . — —— ^ ’ 

*o * 2n w iCn-0 nC«^») 

d 4 = 'A-b/s? , d iU+0 = d tn - Un+l)fe ~^ , 

C 4 - Kl 4® , C zt,H»~ C ai\/Cin+OCZM-4) v 

a(zn-\) 

= l ? ^zn-i — C2-rH»X z<v * 3 ’} 


C. Range of Validity 

Electron self interactions and electron-ion interactions have been 

ne.glected in the present calculations; This is valid for a weakly ionized 

gas, and the purpose of the present calculation is to examine the maximum 

degree of ionization allowable. To do this relaxation times are considered 

to obtain a first estimate of the limit. The ratio of the electron-ion 

relaxation time t . and the electron-electron relaxation time t for singly 
ei ee & / 


charged ions is approximately: t /t . «■ 1/2 /2" 

6© 61 


[46] 


It is sufficient. 


therefore, to compare T ge with x en , the electron-neutral relaxation time. 


The former is given by [47] _ 

zs.g Jv go roif* (fe-T ) ^ 
e 4 n e K 

fa.irC6 0 KT/ g*) 


'-ee 


(III-17) 


where 


A = 


x\ 


V-L 


and r gn is given by [48] : 
Ten 1 


} ^ 4 * 


(III-18) 
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where is the momentum transfer cross-section, U g is the speed of election. 


is tne neutral atom density and E, the mass ratio 2m/M. 


Ilf 4 r» rr f hoco 


results, we obtain figure 24 which shows the fractional ionization when the 

two relaxation times are equal. For this calculation, the momentum transfer 

r 45i 

cross sections were taken from O'Malley 1 , and since £n A is insensitive 

f47l 

to electron density or energy, approximate values for it from' Spitzer L J 

were used. Holt and Hasskel^^ have used a similar curve to define a 

-4 

weakly ionized plasma. They show a critical va,iue of 10 for the fractional 
ionization of helium at one torr and 300°K. There, they considered only the 
directional relaxation. When the mass ratio is added to account for the 
energy relaxation, figure (24) is found to be in good agreement. 

It is noted from figure 24 that as the energy increases, the fractional 
ionization limit is not quite so stringent. At low energies, the charged- 
particle self-interaction tends to enhance thermalization into a Maxwellian 
form. Therefore, if the distribution obtained without the inclusion of 
charged-particle self-interaction is already close to a Maxwellian form in 
an energy region below a few kT, the fractional ionization restriction need 
not be extended to this region. 

For electrons above about 0.5 eV, comparison of the x ge = T gn curve 
with figure (G-4) of Appendix G shows that with a tube radius of 1.27 cm 


and gas pressure of one torr, charged-particle self-interaction becomes 

13 -2 -1 

important for a neutron flux above ~ 5 x 10 cm sec . For gas pressures 
higher than 20 torr, the fractional ionization limit decreases and the 
neutron flux limit can be raised. 

. I 

D. Numerical Results and Discussion 

The choice of some of the parameters required for evaluation of 


Fraction of Ionization 
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Figure 24. Fraction of Ionization Defining a Weakly Ionized Plasma 
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equation (III-16) is not definite, partly due to uncertainties in experi 


mental observations 


r U ^ T - - 1 1 «... 

• rui licxxum etc xvjw cucigito, ojto a 


4- V»o voc*r\Tnl>i 




r>r» r* no 


•F-F-: _ 


-9 3 -1 f 49] 

cient is taken as 4 x 10 cm -sec 1 , the ion diffusion coefficient 

2 

840 cm -tor: 

•16 2 [45] 

om L J 


as D* p = 840 cm^-torr/sec and the elastic scattering cross-section as 


17.6 x 10 cm^ L ' tJJ . With these parameters, Figures 25 and 26 show the 
normalized electron distribution for three pressures. The normalization 
is such that; j f(y) dy = i, where y = . The ion density for a 

O 

given neutron flux and pressure condition is first calculated through the 
procedure of Appendix G. The resultant leakage and recombination rates 
were then used as parametric inputs to equation (III-16). 

A Maxwellian distribution is included for reference. It is observed 
that the distributions shift to a higher energy tail at lower pressures. 

This is because that at low pressures, electrons' are more likely to leak out 
of the system before they can reach lower energies. At higher pressures, 
the energy deposition rate of the heavy-charged particles increases, resulting 
in higher ion densities. However, the' neutral atom density also increases 
causing more electron-neutral scattering which reduces leakage. As a result 
the distribution approaches a Maxwellian form. ' The shape of the distribution, 
not much deviated from a Maxwellian, has also been confirmed by the Monte 


Carlo simulation of B. Wang and G. Miley 


[15] 


Figures (27) and (28) show cases where the tube size has been increased 
2 

to a radius doubling A . Then, for the same neutron flux, the distributions 
for the three different pressures lie much closer together. This is attri- 
buted to the fact that the energy deposition rate by heavy-charged particles 
is essentially independent of the tube size but leakage has been greatly 
reduced by the larger tube size. 

Figure (29) shows an extreme case of changing the tube size at a low 
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Figure 28. Electron Energy Distributions for Pressures of 0.4, 1.6 
and 10.4 Torr and a Tube Radiu. of 1.83 cm. 
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pressure of 0.2 torr. In this case the shift to a higher energy is signi- 
ficant when the tube radius is cut in half. Again, this shift is due to 

the failure of electrons to reach lower energies as leakage increases. 

14 -2 -1 

A more extreme case of neutron flux of 5 x 10 cm sec is shown 
in Figure (30). In this case, charged-particle self-interaction neglected 
in the present treatment, should be included. However, this figure is 
included to illustrate the dramatic shift in the distribution for a high 
heavy -charged particle energy -deposition rate. The peaks shift to about 
1.4 to 2 kT. This is more pronounced in the case of low pressure of 0.2 torr 
where the leakage rate is larger. Reducing the pressure shifts the distri- 
bution in a manner similar to that which occurs as the tube diameter is 
reduced. However, the distribution is not a unique function of dp (diameter 
pressure). This is because, although the leakage term has this form, the 
other terms such as recombinations do not have such a dependency. 

Figure (31) shows the joining of the low energy portion of the distri- 
bution with the high-energy portion developed in Chapter II. The flux 
distribution in the higher energy region is an inverse function of pressure, 
as can be observed from equation (11-10) or from the stopping-power approxi- 
mation of Appendix A. For the energy range 5 kT to 20 eV, a recombination 
r5l) -3/2 

cross section 1 J of E has been used. The collision density, £(E)$(E) 

-3 -1 

cm sec , is a continuous function, but abrupt changes in cross section 
can result in a discontinuity of $(E). Such an effect is observed at 20 eV 
due to the abrupt charge in cross sections used in the two energy regions. 

[col 

Similar discontinuities also appear in other studies and they can be 
eliminated entirely by use of a finer energy division in the numerical 
computation or by adopting smoothed cross sections. 
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Figure 30. Comparison of Electron Energy Distributions for Pressures of 
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CHAPTER IV 


A. The Model 

The total energy distribution of electrons due to an external source 
superimposed on a plasma with an applied direct current electri field will 
be considered. A uniform field, E, is assumed and spatial effects due to 
finite geometry and leakage are not included. The uniform field is often 
good in the positive column of a discharge tube, but it will not be appli- 
cable in other regions such as the cathode fall or the sheaths near the 
electrodes while the field is changing rapidly. Nor will it be applicable 
to some electrode configurations where the field strength varies rapdily 
throughout the system, like the case where a center wire of a tube is used 
as one of the electrodes. 

In moderate to low E/p cases, which shall be assumed, the electrons 
above the ionization energy of the gas atoms are not influenced by the field 
as they slow down. In particular. Appendix I demonstrates that for E/p less 
than 10 V-cm *-mm , the energy gained from the field can be neglected com- 
paring with the energy loss via inelastic collisions. This is the largest 
field considered here and is a reasonable upper limit for other reasons also. 
Most of the applications of interest (see Chapter I) involve smaller fields. 
Also the assumption of a weakly ionized plasma places a similar limitation 
on the field strength, since higher E/p values, like larger external source 
strengths, would result in additional ionization and too large a fractional 
ionization. Finally, leakage is neglected here, and as shown in previous 
chapters, this is consistent with either large systems or high pressures, 
and the latter implies smaller E/p values. 
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With this restriction on E/p, the electron distributions developed in 
Chapter II, for high energy or slowing down region, is therefore applicable 
here. In the lower energy region, i.e. below excitation potential of the 
gas atoms, the energy gained from the field becomes comparable to the electron 
energy. Thus a new treatment of the distribution is required for this region. 
In doing this the electron source due to slowing down from the high energy 
portion and the sink due to recombination will be included. 

B. Balance Equation for the Low Energy Region 

The steady-state Boltzmann Equation with a d.c. elastic field corres- 
ponding to Equation (II-l) is: 




CIV-1) 


r53i 

Following Allis 1 , the distribution function is expanded in spherical 
harmonics. Retaining the first two terms of the expansion, we obtain: 


= {<>(?, cn + ?,(T,UvS/u CIV-2) 

where, U ) = 

(t,u) =r J-fCT,u)udrL 

Substituting equation (IV-2) into equation (IV-1), one arrives at two 

rs 3 1 

coupled equations : 


± - et- 

r^io.cf.u)\ = fs^ 0 \ + rs£ tt \ 

' $"t_ JG- \ b~t- J eiajji it S"t_ l inelaiiic . 


(IV-3) 

(IV-4) 


''■*1 JeUshc u V t lTo r* e ° 


where 
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= the net rate of influx of electrons due to elastic collisions in 
the velocity element d^u 

vSt /inelastic = the net rate of production (loss) due to inelastic 


3_w 

collisions in the velocity element d u 


(f| ) = rate of momentum influx into the velocity element d\i 

T = gas temperature 

and v c = momentum transfer collision frequency. 

If the first-fundamental-mode diffusion approximation and uniform spatial 
properties are assumed, f(r,u) becomes spatially independent. Then, combining 
equations (IV-3) and (IV-4), we obtain a second order differential equation 


for f (v) : 
o 


S ■+ 


(A 


U 




\-C + -L- ± G - O 
4 uu' C*u ^ 


(IV-5) 


where S = total source rate from primary 6 -electrons and ionized secon- 

3_ 

daries into the velocity element d u. 
y r = recombination frequency 

A = characteristic length of the system from diffusion approximation 
and G = — 4JT ft, -J- d_ C _x 3. m \ 

3 U v c hu gyV, ■+ ) 


with y = /'eE'f'J- 
c ^ 

Omitting the source and neglecting spatial losses and recombination 
while assuming [<*'/(m\ + M)] k~f , Nigham^ solved equation 

(5) numerically for CO and CO 2 gases. However, for present purposes, a 
formal expression for f can be obtained from equation (IV-5), (see Appendix J) , 
as : 


U 




A(u‘) 


J 


(IV-6) 
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i 

“here BCU ^ Km\Ac. 


and 


A CO = 

S'(u') = 

X(c) = 
a(c) 


F z c AcO + mkTc . 3 

3 M\Cc) 


• u 


^ [r) a £ °CO -c< CO*] O Ac. 

J o 

c A)CO ^ mean free path 
electron production rate in velocity element d~u 

m 


3- 


The assumptions of (m+M) = M and no leakage losses have been made. 

Equation (IV-6) can be solved by a simple numerical integration provided 
that an approximate form of f°(u) can be obtained for the calculation of 
S'(u'). The case without external electron sources or recombination has 
been solved by Druyvestyen ^ ^ , Smit ^ ^ and others Their results 

have 'been used as the initial estimate of f°(u), and a corrected form of 
f°(u) is obtained for the present case of interest where a source and sink are 
present . 

This procedure was used in preference to the direct numerical solution 

of equation (IV-5) for two reasons. First, a Monte Carlo simulation for the 

case with an electric field shows that the distribution is rather close to 
T541 

that of Smit 1 . Second, the numerical solution of equation (IV-5), such 
as used by Nigham, is also a first order approximation since a two-term expan- 
sion of f(r,u) was used but eventually only the first term was calculated. 


C. Numerical Results and Discussion 


The uncorrected f(u) for no source or recombination, similar to that 

r 2 1 

discussed by Deuyvesteyn 1 J ,was used to estimate the rate of recombination 
and rate of ionization. With the assumption that 1 ^ , 


64 


f(u) becomes a function of \CO or E/p. Following the work of 
. . [54] 

smit - , ' , a(cj is assumed to be a constant for electron energies below 5 eV 
while from 5 to 19.8 eV, A(c) is assumed to be proportional to c. From 19.8 
to 24.7 eV, Smit has obtained f(u) represented in Bessel functions. 

r 2 1 

Druywesteyn et al 1 1 show that these Bessel functions can be approximated by 
an exponential function and this simpler form for f(u) has been used from 19.8 
to 24.7 eV. 

the total electron density n , is first estimated by substituting the uncor 
rected form of the distribution into the appropiate rate equation, namely: 

CIV-7) 

-u 

where 


ao 

-p/w - fte ( f( E/P.t) <T; Ce) At = O 


B = recombination coefficient 
p = rate of energy deposition by heavy-charged particles 
W = energy required to create an ion pair 
f(E/p,e)de = normalized distribution in energy interval dE at e 
cr^(e) = ionization cross section from Vriens. 

The first term of equation (IV-7) represents the total rate of recombination 
sink S^, and the last two terms represent the total rate of electron pro- 
duction S. To compute the differential source rate a(u)du, we use the 


differential energy loss cross-section o(AE) of Stabler 

4TtACU)U 7 'dU = ■= 5 crCx)c\x 


[25] 


such that : 


(IV-8) 


where x = AE-u, and the normalization is such that : 

\ S (TCX) AX — s . (IV-9) 

■'o 

For computation of the differential recombination rate, the uncorrected 
form of f (u) is again used. Equation (IV-6) has been evaluated with these 
sink and source rates to obtain a first order correction to the distribution. 
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Numerical values for the source were selected to represent typical ionization 

rates due to heavy-charged particle energy deposition. Figure (32) illustrates 

13 -1 -1 

the distributions for a neutron flux of 5 x 10' cm sec in E/p = .3, 4, 6 and 

-1 -1 (* S 

8 V-cm -mm . The nomalization is such that \n(E)dE = 1. The population over 

J o 

25 eV is considered to be negligible in this normalization. Figure (33) illus- 

14 

trates the distributions for the same E/p values with a neutron flux of 5 x 10 
- 2-1 

cm sec . The deviation from an uncorrected Druyvesteyn-type distribution is 
largest when the neutron flux is high and E/p is low. Two such cases are compar- 
ed with the unperturbed distributions in Figure (34). In general, the perturba- 
tions are not too large and the use of unperturbed distributions in the source 
and sink integrals is justifiable. . This is especially true in the case of the 
smaller neutron flux with E/p larger than 4 V-cm *-mm 

The perturbation is examined more closely in Figure (34) which compares 

the E/p = 3 and 4 V-cm ^-mm * cases for no source or sink with the cases where 

14 -2 -1 

the neutron flux is 5 x 10 cm sec . It is seen that as E/p increases, the 
deviation decreases. This is because in higher electric fields, the distribu-r 
tion is more dictated by the field than the electrons slowing down from higher 
energies. The ratio of electric field strength and the source strength, th^n, 
determines the importance of the inclusion of source and sink in the calculation. 

The average velocity u can be computed from the distribution obtained. 
This random velocity is also proportional to the drift velocity . The 
average electron flux u n g) therefore is proportional to the current density. 

A plot of the average flux vs. E/p is shown in Figure (35).. This has the same 

[571 

characteristics as the V-I measurements of Ganley 1 ' J in the enhancement of CO„ 

L. 

lasers by nuclear radiation experiment. Although an exact comparison cannot be 
made due to the different gases involved, the trends in Figure (35) do agree 
with the measurement in CC^. 
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Figure 32. .Normalized Distributions in Electric Fields , E/p = 5, 4, 6 and 

8 Volt -ca ^and a Neutron Flux of 5 ?. “-sec * 
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13 -3 -1 

It is concluded that for neutron flux of less than 5 x 10 cm -sec 

i 

(_• equivalent to a primary source of about 5 x 10 1-keV electrons in 

10-torr helium) the distribution is close to that arrived at by Druyvestyen 

or Smit. This has also been demonstrated by a Monte Carlo calculation ^ ^ 

where only moderate deviations are observed only for primary electron (1 keV) 

16 3 1 

source rates of 10 cm ,-sec or above. It should be noted, however, that 
in the case of helium, excitation and ionization potentials are relatively 
high. Larger perturbations can be expected in the cases such as molecular 
gases, where lower energy levels are involved. For those cases where devia- 
tions become severe, equation (IV-6) is not accurate unless a better estimate 
of f 0 is available for the source and sink integrals. A numerical solution 
of equation (IV-5) would then be necessary. 
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CHAPTER V 

CONCLUSIONS AND SUGGESTIONS FOR FUTURE EXTENSIONS 

A . Summar y of Present Work 

1 . General Remarks 

The energy distribution of electrons in helium resulting from a 
continuous high energy source has been analytically studied. The main 
purpose of the study is to provide a method to analyse the problem. Within 
the illustrated constrains, such as the degree of ionization, and with 
appropriate changes in reaction rates, the method can be extended to cover 
other gases or gas mixtures. A set of physical parameters, such as dimen- 
sions of the system, gas pressures, and primary source rates were used to 

illustrate numerical results of the analytic methods developed. These para- 

[13] 

meters were consistent with experiments of J. Guyot et al. in studying noble 
gas plasmas created by nuclear irradiations. Again, within the specified 
constrains, these parameters can be changed to suit new situations as they 
arise. 

If analytic results are not feasible when new parameters are intro- 
duced, the basic equations of the present study are still valid and may be 
tackled on a purely numerical basis. As shall be detailed later, the present 
work complements the Monte Carlo simulation of the problem and it may be of 
advantage to incorporate both methods in a single program. However, some 
of the specific results of the present analysis shall be pointed out first. 

2 . The Analysis of the Distribution in the High Energy Region 

For energies above the excitation potential of helium, the inelastic 
collisions were found to bo the dominate process in the slowing down of elec- 
trons. The balance equation has been set up for this region with the approp- 
riate slowing-past kernels calculated (equation 11-28) . These kernels were 
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derived from energy-loss cross-sections. The solutions for the flux dis- 
tributions were numerically obtained from the balance equation. They are 
critically dependent on the slowing-past kernels and therefore accurate 
cross-sections are necessary. The slowing-down past kernels and appropiate 
limits developed here are uniquely suitable for this energy range accounting 
-for the secondary electrons. 

The W-values in helium were calculated to check the validity of the 
flux calculation and to determine the appropiate cross-sections to use. 
Vriens' energy-loss cross-sections were used. Emperical parameters were 
introduced to force the integral of these cross-sections to agree with ex- 
perimental total ionization cross-sections. The W-value calculations alone 
do not clearly demonstrate the superiority of the fitted cross-sections but 
they do indicate that when individual excitations were considered, the agree- 
ment is better than the case when a one-excitation level model was used. W- 
values calcuated for a 1-KeV primary source using the six largest excitation 
cross-sections were 41 eV per ion pair using the fitted Vriens cross-sections 
This is to be compared with the W-value of 42.3 eV per ion pair measured by 
Jesse et . al . 

The effect of leakage was included through a first-fundamental model, 
diffusion-approximation - leakage term. The effect of change of pressure and 
system size is observed, and the shift in the distribution also clearly 
illustrates the effect of leakage (Figures 21-25). It is shown that the flux 
depression due to leakage is more severe for higher energy primaries, for 
smaller systems and for lower pressures. . 

These flux distributions cannot be obtained if only the therr.alized 
swarm is considered. For a source rate of iO^cm "^-sec ^ , 1-KeV primaries, 
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the Maxwellian tail at 50 eV would underestimate the flux density presently 

-4 

calculated at that energy by a factor below 10 . If the approximation of 

dE 

flux ~ const Jl Tj- > were used, it would overestimate the flux density at 
50 eV by a factor of more than 7. The flux distribution presently calculated 
would therefore provide a much more accurate electron reaction rate than 
other crude approximations. 

3 . The Distributions at the Lower Energy Region 

For low energies, a model similar to the "heavy gas" model in neutron 

physics was found to be applicable. This model is uniquely suitable for the 

presence of a high energy electron source. The fraction of ionization or the 

primary source rate beyond which charged-particle coulomb interaction“cannot 

be ignored has been determined through an analysis of relaxation times. For 

a 1.27-cm radius boron-coated tube at 1 torr, this limit corresponds to a 

13 -2 -1 

neutron flux of 5 x 10 cm -sec 

For the cases without an electric field, the distribution is found 

to be close to a Maxwellian with deviations due to leakage and recombination. 

A shift to a higher average energy due to the inability of higher energy 

electrons to thermalize is noted. For the cases with an externally imposed 

direct current electric field, it is found that the high-energy electrons 

are not appreciably influenced by a moderate field of E/p 4 10. The low 

energy distributions are close to the results of Druyvestyen when no external 

sources are present. However, for a large primary electron source rate, 

14 -2 -1 

corresponding to a neutron flux of 5 x 10 cm -sec , deviations from the 
Druyvestyen results can be observed. Finally, the enhancement of current 
densities in the presence of a primary electron source are noted. 

For cases both with and without the field, the model results in a 
second order differential equation. Solutions were obtained by a series 
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expansion and a perturbation approximation. Numerical solution of these 
equations may be more appropriate in cases (for instance, in the presence 
of molecular gases) where the algebra becomes too involved or when the 
deviation from a Maxwellian or Druyvestyen distribution becomes too large. 

Many of the results obtained in the present study agree with a Monte 
Cprlo simulation reported by B. Wang, and these are pointed out more specifically 
below. 

4 . C omparison With Monte Carlo Results 

[ 32 ) 

The Monte Carlo simulation 1 J and the present work show the same 
general shapes in the high energy region. In the presence of an electric 
field, the Monte Carlo result show that for the case of E/p = 10 V-cm ^-mm ^ , 
the distribution in the high energy region is not influenced by the field. 

For the low energy region, the distributions can be approximated by a 
Maxwellian or for the cases with fields, by a Druyvestyen distribution. 

Deviations from those were observed when the source rates are extremely high. 
(>10 A cm -sec , 1-keV electrons). These deviations have the same trend 
as the deviations cf the present calculations relative to source rates, pressures 

5 . Concluding Remarks on the Methods 

Numerical results under various specific parametric conditions were 
calculated from the equations developed for the entire energy range below 
the source energy. W-value comparisons and comparisons with a Monte Carlo 
simulation show these specific results to be in good agreement. Comparisons 
with temperature measurements by micro-wave techniques I" 58, 59 ] are nQt mac je 
because in the presence of a high energy source, the thermal ized condition 
is not assumed. As indicated, the important high energy region is indeed 
found to be highly non-Maxwell ian. The agreement of current enhancement in 
the presence of a high energy source with Ganley’s^ ^ observations has been 
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noted although no numerical comparisons were made because of the differences 
in gases and reactions involved. 

The basic equations are generally valid for other gases and physical 
parameters; however, the accuracy of results depends heavily on the cross- 
sections and kernels. For the high-energy region, equation (11-28) can be 
used provided slowing-past kernels are derived accordingly from accurate 
cross-sections. For the low-energy region, equations (IV-5) and (III-ll) 
account for the presence of source and sink for the cases with and without 
an electric field superimposed respectively. These equations can be numerical- 
ly solved without any of the approximations used in the present analysis, i.e. 
the series expansion used in solving equation (III-ll) and the perturbation 
method used in solving equation (IV-5). In fact, these techniques were success 
ful in the present case partly due to the specific circumstances involved, such 
as cross-sections in helium, pressure ranges, electric field strengths and 
system sizes . 

B. Future Extensions 

It is suggested that future work should concentrate on applying the 
methods developed here nad in the direct confirmation of their validity 
through experimental observations. Some examples are: 

1. The electron energy distributions can be calculated by the methods 
developed for the appropiate gas or gas mixtures. These distributions can 
then be used to find electron reaction rates for a set of rate equations 
describing the plasma kinetics. The emission spectrum from a certain excited 
state, for instance, could then be measured to varify the calculated excited 
state density. This provides an accurate check of the calculated electron 
energy distribution. 
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2. The energy distributions also provide a starting point in several 
applied calculations in a radiation-induced lase'r application. The electron 
energy distributions can be used to determine possible inversion of excited 
states by solving the rate equations described above. In a nuclear light- 
bulb reactor, the appropriate distribution can be used to calculate the 
radiative emission spectrum from the fuel region and from the buffer-gas 
region. This sets up a criterion for the addition of seeded gases to prevent 
harmful radiations from reaching the wall structure. However, for a compli- 
cated situation such as the gaseous core reactor, Monte Carlo calculations 
may be a valuable supplement to the present methods. 

3. The present method can be used to complement a Monte Carlo calculation. 
Because of the small electron population in the high-energy region, in a Monte 
Carlo simulation, a proper treatment of this region results in a relatively 
large number of sample particles in the lower energy region. Weighting 
factors may help, but the calculation may still be prohibitively expensive. 

The present analytic method handles the high-energy region with good efficiency, 
but involves complicated expansions and approximations in the low-energy 
region. Thus there is a strong motivation to combine the methods. The 
present method would be used to generate a distribution in the high energy 
region while a Monte Carlo calculation would provide more accurate results 
in the lower energy region. 

The present calculational technique might also aid the Monte Carlo 
calculation by providing an initial input. The convergence of a Monte Carlo 
code depends on an estimated initial input. A good initial input can be 
obtained through the present calculation to assure a rapid convergence of 


the code. 
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APPENDIX A 


Approximate Flux Using a Continuous Slowing Down Model 

With No Secondaries 


U. Fano-^', using transform techniques on the flux conservation 

equation, demonstrated that to a first approximation, the flux is related 
* 

to the inverse of the stopping power. This can also be shown by expanding 
the flux in the particle balance equation. In the absence of secondaries, 
the balance equation becomes; 

SCe'JICE') = j + S^Ce-Eo') (A1) 

0 

where 2LC&. 1 ) is the total collisional loss cross-section, kCE +1, I) 
is the probability of a particle of energy (E'-i-t) losing an energy of 
amount t from energy, and SpSCE-Fo), the primary electron source, is taken 
as a delta function at E q with source strength Sp . 

If $(e) and k(E,x) are slowly varying functions of energy, the first 
two terms of a Taylor expansion give: 

?Ce 1 ^Ce.") = + ^ Ce ’^k(E.',OT:clx+S f SCe-E0 (A2) 


In the absence of absorption, 

zl 


\ kU',T)dT. -=• kU'x') * H(e') 

O Jo 

Multiplication of equation A2 by dE* and integration yields: 

E -CO 

° = fldt dE’ i-.^CE’jkCEU) + s f 


(A3) 


and since 


inn $CE)k.CE,'E) = ° 

E -*«o 


(A4) 


(A5) 


we obtain: 


$(E) - c.» t/ (' We . >7 . ix 


(A6) 


[4 1 

Equation (A-6) agrees with the result U. Fano L 1 obtained by Laplace 
transforming equation (A-l). 
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APPENDIX B 


Derivation of an Integral Form of the Collision Density Balance Equation 

After an ionization collision, two identical electrons will emerge. To 
have an accounting system we shall label the one with larger energy as the 
"primary" electron and the one with the smaller energy of the two, the 
"secondary". The collision density balance equation for a 6-function source 
of electrons at energy E is : 

H c (&) $ (e') = \ -V S^ + S SCe-e*} (B-1) 

E-+ u, 

where 3 Z^e'-*E ) is the differential cross-section of scattering 

into energy E from energy E'. 

U, = first excitation potential of the gas atoms 

A =E+Uif2E+U$E , and A = E - E if 2E + U > E . 

o’ o o 

The "secondary" source term takes the form: 

S s Ul = ^ X s ( £.'-•> e ) ^Ce^cIe 1 C b -2) 

lE + O 


where X s (e‘ t ^ ie probability per electron track length that the 

"secondary" electron emerges with an energy in dE at E. 

We shall divide the energy range into two cases: . 

Case 1. If 2E + U E , then S (E) = 0 and A = E - E. 

o’ s v 1 o 


Integrating equation (B-1) from E to E q , we obtain, 

\ - ( E ' dE "f° H c CE'-.E-)$CE , )dE' + S p 

e i P 

Now consider the double integral on the R.H.S. of equation (B-5) . Since 
for E" in the range of E q - to E^ , (E 1 — > E") = 0, we have 


(B-5) 
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V- — ( \ t 

“ r 


X c Ce!-* e )c!e' 

e"+o ( 

(B-4) 


*= ( dt" 

E J EVU, J e 

r E ' + u, 

= c 1 e‘ $(E‘) \ dE" X«X E 1 -? 

E-»U, J E 

where we have interchanged the order of integration. Combining equations 
(B-3) and (B-4) , we obtain 

fE'-U, pt. 

j ^CEUE^dE, = J dE'^CEoW.ICfUE.l+Sp ( B ' 5 ) 

e EtO, ' 

where E . is the minimum energy above which electrons with energy E' can 
nun oJ 

slow down. Note that we have changed the notation of dummy variables 
inside the integrals. Breaking up the LHS of equation (B-5) gives: 

/-e-HA pE'-u, r E » pE'-f, 

LHS = j $C^')dE'J X^E'-^E^Ed JdE'fCE'H Z c &e)Ae x (B ~ 6) 
e Ern "* E+U, F miv) 

so equation (B-5) becomes 

E+O, z SO, r E " U ' 

1 $CE l )de' C + \ He'^Ce 1 ^] ’Z c CE , -»OdE, 

J p Jc ^E+L) E*nnJ 

' - ^ X c CE , -?E t )dE l J - S f 


OY 


(B-7) 


r ^u t rE'-U, r f 

J <£(&')c\fc'j X c (EUf,)clE ( -v 1 $ (e')cIe 1 1 ZcCe.Lt, )oIe, s ^ 

Emw r 1 , * . P 

c+U, c mir. 

Consider the L.H.S. of equation (B-7). The flux 4>(E') in the first integral 
lies in the range from E to E + U, and any inelastic collision will slow the 
electron down past energy E; as for the remaining integrals, E^ is in energy 
range E . to E and E (F. 1 -*■ E.) accounts for all collisions that slow the 
electron down past energy E. Therefore, the LHS of equation (B-7) represents 
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the total rate of slowing down past energy E. 

Case 2 . If (2E + U) < E q ; in this case A = E + U and 

r e ‘ 

S s c^=\ Y^CE^E") 

^ze-*o 

Integration of equation (B-l) from E to E q gives: 


(B-S) 


E 0 Eo-U -lE VO 

( “l c (E'^ CE"^E" = ( 1 at*' dE' §Ce')ZJe'-=>E") 

J E J E Je'Vo, 

+ j J j dE'f(E>)I c CE'^EV H J S s C^c|E‘ -v S p 


(B-9) 


Consider the RHS of (B-9). The first two integrals cover the areas (1) 
and (2) in the sketch of Fig. B-l. This is equal to areas [ (1) + (2) 

+ (3) ] - (3), i.e. the first two integrals become: 
e 0 rE'-u, - e„ - eU/ 


r So r t-w, ^ 

J SU'ldtM dE'^^EVE 11 )- <§(E^dE'\ cA<="^Ce'^EL") 

E+U - E ie-t-U E. 


(B-10) 


But, for the second integral, Z c (E' -»■ E") is equal to zero for the region 
of E" since no "primary" electron can slow down past (E' - U )/2. The 
integral therefore vanishes in area (5). Then equation (B-9) becomes: 


V f E. r 6 -U, E 0 

\ Z t CE‘ , ')$(E' , )(iE ,, = J $CE')dE'\ 4 e u I c CeVe'') + j 5 /eW -w 

l e + u. e e 


(B-l 1) 


Equation (B- 11) is identical to equation (B-5) except for the "secondary" 
production rate term. Using equation (B-7), equation (B-ll) can be 
written as : 

r E + U, ,E'-U, rE 0 rt= 

) ZcCEUEOdE, t \ ^CeVdE' \ ^CeU^dE, — 

6 *** J ErU, L 


+ f S s C&')dlE' 


(B- 12) 
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Here the LHS is equal to the total rate of slowing down past energy E while 
the RHS equals the total rate of production above energy E. This is the 
same as: 

- r E “ r Eo 

X- \ K-CEj e) dE' — ^ \ s (EOdE 1 (B-13) 

*1 E J E 

where K^(E',E) = probability per track length of slowing down past energy E 
through a j th type collision. 

This shows that in both cases the differential form of the balance equation 
(Equation (B-l)) and the integral form of the balance equation are equiv- 


alent . 
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APPENDIX C 

Comparison of Average Energy Loss Due to Elastic and Inelastic Collisions 

For High Energy Electrons in He 

For elastic collisions, we consider the first Born approximation 
[59,60] 

cross-sections . 


r 00 

IC©'] = | 1 JLjft i* j 4in k-t V 


CC-1) 


where k = .2.M c O .■ p./ 

— Am 0 /*. 

K 

For helium, Hylleraas^^ has employed a variational method to obtain the 
following wave function: 

, , . x -iCr.-k-Y^/a 

m -OVtro.Oe tc-2) 


with the potential: 

V O') = 


- ie_l (A + 1A e 


-zEr/<x< 


(C-3) 


where z = 1.69, and a is the Bohr radius. From this, the differential 

0 * 

scattering cross-section becomes 


T(f o, = CjjtzAlX 

1 


(C-4) 


where A = 2 z/a , A = atomic weight. 

• 0 

The fractional energy loss due to elastic collisions is related to 

1 


the angle of scattering by 


4! = _ cos ©) 

£ M 


(C-5) 


where m/M is the mass ratio between the colliding particle. From equations 
(C-4) and (C-5) the average fractional loss per track length becomes: 
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\ 4 nA 1 K 
) \ + 

f ' 

S x 3 

** 0 

Oo^+ 

(a. -t- b x 1 -)^ 

- dx 


- S" Fpr) 

4-ri^e^ A. 

r - 1 - j 

? a + t 
■3 b + 

1 

a 

It 4 - 

_ v- ’ 

4b" 

bK*+b) 


o> 


4b’ J "C0-+b)' 1 bkHdt-b) 


-°£ 1 

) x Cd + b ) 3 J 


where a = \ , 

b = c-^X = 

After rearranging the constants, we obtain: 


j 1 = L»*X 

iT ) v- 


<*1> - {hl' + (>T 3 ; -k - -F 




(C-6) 


*v 4- 

^ . + 

k 1 


4 


bC^+ySM 3 i 

(C-7) 

u /c. 

and It = ^ C 

V*iT 



and we have used the fact that X’ « -t 

-w. 

For inelastic collisions, we consider the Moller cross-section 


i rm ^ ^ 


L(E ( ^dl = rn]V^ f X + J- - ^xkL[ 1 v _L_ 4 - J— 

p II V (F'^ <e + 0 L l CE+, ') al I(C-8J 

where E and z are the incident energy and energy loss for the collision 

? 

respectively, in units of M^C - . The average fractional loss per track length 


becomes : 


<v>,„ - j 

l 

Q 


\d-L 






\ 1 ^ T 

(. 7 -+ E ) 



k 


l X (E-t) 1 - 

(be] 1 


% 

[ 

1 

X 1 

+ X l 





E--L J 

0+ E)- J 


(C-9) 

2 2 

where r = e /me", 
0 

and Q 

is a constant determined 

from stopping 

power 


theory. For U << C, Q equals Bethe's^J ^min = ^ | Z ' ~o • The upper limit 


is arrived at from the consideration of indistinguishability of electrons. 
Changing the energy units to electron volts and integrating (C-9) , we find 



3.H Me 6* 
E u 



& 

E - Q 


E 

ZQ 





E/t 

E-Q 


+ j . . E A- ^ _J f|’_ s l )l (C-10) 

(.E+mc'-r J E-Q Ce-hwc'^ 2- 

Comparison of equations (C-7) and (C-10) demonstrates that in the energy 
range where Born approximation and Holler's cross-section are valid, the 
average energy losses due to inelastic and elastic scattering differ by a 
factor of magnitude about m/M. This justifies the neglect of elastic 
collisional losses in that energy range. 
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• nnr\:mv r\ 

ftrrcftuiA u 

Slowing Down Kernels from Binary Encounter Collision Theory 


The binary collision theory neglects the influence of the atomic 

nucleus and treats inelastic collisions as an interaction between two 

electrons alone. IVhen two interacting particles obey the inverse square 

law, the problem is exactly soluble in wave mechanics, and the solution 

r 39 1 

yields the same scattering law as the classical theory 1 J . If, however, 

the two interacting particles are identical, as in the case for electrons, 

symmetrical and anti symmetrical wave functions in the space-coordinates of 

the electrons must be used. Therefore, if the corresponding cross sections 

are <T + and CT and ® is the angle between the spin directions of the electrons, 

[391 

the total cross section is L 


~ 4 ( I - Cera 0)0"^ ^ C 3 t un © ) O 

If one assumes unpolarized electrons, this reduces to 


CD-I) 




(D-2) 


L. Vriens' 
interference gives 


model of tw ; o interacting electrons with exchange and 

[ 20 ], 




ire 4 



4- 

3 



_J 






CD-3) 
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^4 r~ / \ 

^ d4E - X u 


A C. \ 


\ \ + 
4E 3 


4E* + | CE*-E X - 


4E) 2 


+- 


4 E4 


^le v E a -AE\- 




e! 


AeiEi-t^-AEl 


-1 ( w 


, ?-«• El $■ (D-4) 


where the terms are defined in the text. The value of $' is ~ 1 for 
E 1 >> R. In the following integrations, we shall approximate $' by a 
constant. This is consistent with the assumption that Vriens made in 
arriving at his total ionization and excitation cross-sections. 

For the case of ionization collisions, consider E^ to be the incident 
energy and E^ - AE and ZiE - U to be the exit energies after ionization, 

U being the ionization potential (Fig. D-l) 


E, 



e.-^e 

^AE-U 


Fig. (D-l) Electron Energies in 
an Ionization Colli- 
sion 


The two emerging electrons are indistinguishable and the maximum energy 
loss possible is' (E^ + U)/2. The minimum energy loss that an incident elec- 
tron can suffer and still slow down past energy E from energy E' is E' - E. 
The siowing-past kernel due to ionization collisions is then obtained by 
integrating Vriens' differential cross-section through these limits. Thus, 
for (E' + U)/2 >_ E' - E >_ U; 

AE'+'-'J/z. 

Me' ( E) * J dAE 

e‘-E 


v4u^ 4 

*> - 


» 

Ef'-E 


__l 




L_ 

-Ce'-eT 


-+ 


i 


] 


(e'+ u) 
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where 


—a 

$ 


(Lets 


r fo~ 


u.-E.y 


u 


■E'-E 


(D-5) 


and N g is the total number of valance electrons in helium. For more compli- 
cated atoms, one has to sum over each different electrons in the shells with 
their appropriate ionization potentials and averaged kinetic energies. An 
approximate method suggested by Vriens is to simply multiply by the number 
of electrons in the outer valance shells. 

If, however, E' - E <_ U, any ionization collision will show the electron 
from energy E' past energy E. Therefore, for (E ' - E) <_ U; 

r ESu 

KxCe',E^ - \ ^U e C5; E ciAE 
u 



This is, as one would expect, identical to the total ionization cross-section 
and is consistent with the ionization cross-section given by Vriens and 
others. The two slow r ing-down kernels (D-5) and (D-6) are continuous at 
E' - E = U. For E' - E >_ (E' + I)/2 or E' <_ U, the slowing down kernel 
Kj (E ' ,E) = 0. 

The binary collision model assumed that an excitation to level n 
takes place when the energy loss AE is such that, U < AE < U ,, where. U 
and are excitation potentials of levels n and n+1 respectively. In the 

case of helium where the final state depends on the spin orientation, we 
w'ill 3pply the following approximation. Mott and Massey give the excitation 
cross-section for an arbitrary angle 6 between spin directions of interacting 
electrons to be^^: 


Q e (© , | ~ Q de +Q ee ^ (»-«•<«» B) Qi«t . 


CD-7) 
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where de, ee and int denote contributions due to direct, exchange, and 

interference. For parallel spins, Q (0) -» Q , + Q + 2 Q. and for 

c de cc mt 

antiparallel spins, Q O) = Q , + Q . However, direct excitation only 

0 Q 6 C6 

occurs in the singlet mode while exchange excitation can be either singlet 
or triplet. If one assumes one half of electrons have opposite spins and 
the other have parallel spins; then 


Q e = \Q e C^) -*• iQ e Co-) = Q^A) 
— Qde + Q ee. Q.'nt 


CD-8) 


From this we deduce that, 


* Qde ^ % e 
In general the excitation to the n^ 
to: 




"taipLei 



(D-9) 


cross-section is calculated according 


Qe.Cn) 



(D-10) 


Since- only AE, the loss of energy is considered, there is no distinction 
between levels when one considers the kernel of slowing down past a certain 
energy. One can therefore group all the excitations into one level and 
consider an excitation range, i.e. one considers excitation to have occurred 
if Uj — ^ — ^i w ^ ere is the first excitation potential energy and IL , the 
ionization potential. We refer to the slowing down past kernel due to exci- 
tation derived through this model as the one level excitation model. 

For E' - E >_ , no excitation collision can bring an electron from 

energy E' down past energy E and the slowing down past kernel K (E',E) = 0. 
For E' > U, and U. > E' - E > U, 


91 


Ks*. 


Ce'.e'i - \ 'N e CT. E i A s= 

e'-e 
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1 J 


L(e'-E) 1 11* E 11 - (.E-tllf. 

IeVo.} 

L (E'-E)E' J 


where 


(D - n) 

For E' >_ and >_ E' - E, any excitation collision will bring the elec- 
tron down past energy E from energy E', and 


U- 

K«(el^ =5 AaE 

U , 


r 

( J + 

f± _ 

' )>'**>■ 

CE'+vaJ) L 

uj + 

l E* 

e‘ ■+ u. - l?, 1 "* t> 

-L + _i_ . 

_ 1 1 

< 

i rucCE'^-^n 

U> 

Ce'+u.- u.y-j 

E'i-o x ; 

"1 E'U V Jj 


CD-123 


This agrees with Vriens' one-level total excitation cross-section. 

The secondary production kernel is defined as the probability per track 

length that an electron with energy E' is going to produce two electrons 

N \ 

with energies both larger than E. This (Fig. D-l) imposes the requirement 
that AE - U > E or E' > 2E + U. For energy E' >_ 2E + U, after an ionization 
collision, at least one of the electrons will emerge with energy larger than 
E. For both to emerge with energies larger than E, the minimum energy loss 
of the oide.nt electron must be E + Li. The maximum . loss of energy due to 

v-.stinguishibility of electrons, cannot exceed y (E 1 + U) . The kernel 


the in 
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for secondary production K (E',E) is; for E' > 2E + U; 


A(e'-t-u') 

K s IE',£)=] 

EtO 

= < _\ ! liv r 1 

(e,+»0 Ie-vO E ‘-& s L Ce + uf- 



(D-13) 

and K s (E',E) = 0 for E' < 2E + U. 



APPENDIX E 


Determination of the Number of Secondary Electrons Produced Per 6-rav 


First, we shall show explicitly that in an energy range below the first 
excitation potential but above an energy where recombination is significant, 
i ; e , between 19.5 and say 2 eV in helium, the slowing down past density q(E) 
is a constant. In this region, the recombination is considered to be negli- 
gibly small and in the case of helium, processes like dissociative ioniz- 
ation (e + A 0 -*■ A + + A + e) can be neglected ^ ^ . Consider the collision 
density F(E) defined as the total number of collisions at energy E per unit 
energy interval; 

r u > 

FCe.') = \ F(e‘)P(t l ^E)aE'ciE -S^dE (E-!) 

. J e 


where 


FCe 1 = Z s Ce1SCE) , 
= , 


and Sj(E)dE gives the rate of electrons arriving in dE at E from inelastic 
scatterings above energy U . E (E) and E [E) are macroscopic elastic and 
recombination cross-sections and PTE’ - E)dE is the probability that an 
electron scattering at energy E 1 will arrive in energy interval dE at E. 
Integration of equation (E-l) from 0 to E where E is below gives 


e f e 

^ FCe^e' 1 = j FCeVtE , ^E")dt' + ] s a ce."i] 4E 1 ’ 


(E-2) 


The shaded area in sketch is the area of integration for the double integral 


Figure E-l 
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on the R.H.S. of equation (E-2) . Changing the order of integration, we 
write the double integral as a sum of integrations over the two shaded 
areas (1) and (2); i.e.; 

f de"f c\E* = 1 FCe 1 ') f pCe'-^E") 

J ) j— ii J „ 


f o, ,* 

+ \ dE'FCE')] PCe'-9E W )c\E‘' (E-3) 


But 


•r 


- i . 


(E-4) 


and 


\ PU'^e.'OdE" ) 


where K(E',E) is the probability that an electron slows down past energy E 
after collision at E'. From equations ( E-2 ), (E-3) and (E-4); 

c E c E 

j F(E“)dE" ~ ^ de'FCE^H- J S x CE)cke' 


-V J Vf^KCSL' E) , 


(E-5) 


or 


U. 


dE'F(E’)^Ce;e) ■+ S T Ce)dE v = 


e , 0 o (E-6) 

But the L.H.S. of equation (E-6) is exactly the definition of q(E); i.e. the 
rate that electrons slow down past energy E or appear at energies below E. 
From this; 

r c 

s~\ r \ l <r *- • 1 — I 

(E-7) 


^ s^ce'^dt 1 


However, above a certain energy E^, recombination is not significant, so 
we obtain 


0 


(E-8) 


. L. 
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and q (E) is not a function of E in this range. This value of q(E) represents 
the total number of electrons slowing down per unit volume and time and can 
be used to calculate the W-value of a gas. From equation (E-6) 

r U| 

5 x C&'\dE' 

"O 

= f U, dE' ( 5Ce")r, M Ce' , -’E'i 

h, 

= r*<\6’$U“) KCe",u,v 

V 

CE-9) 

where o 

J o 

= number of electrons produced with energies 
below per unit 3>(E") at E" around dE". 

Using the previous notation for slowing-past kernels due to excitation and 
ionization, we can illustrate some of the K^ n (E",U^) in different energy 
ranges : 

For E" >_ + 2U^ , after an ionization collision only one electron 

could have energy less than U^. No excitation collision is going to result 
in slowing down an electron past energy U^. In this case: 

e^u.-o.-i (E . 10) 

where = U, the ionization potential. 

Using similar lines of argument, we obtain 

U. + 2U. > E" > U. + U, : 

1 1 — — l 1 


Kt* Ce'.'u,') ^ KjCe", E'-iU + CE . n) 
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U. + U. > E" > 2U • : 
1 1 — — 1 


, (E-12) 


21L > E" > U. : 

1 — — l 


K.kjU’'^= zk x CE",E ,l -O.J + t 

U. > E" > U. : 

1 — — l 


k^cc‘;u,) - 

U, >_ E" : 

kCe". U.) - O 

irsi 1 


(E-13) 


(E - 14) 


(E-15) 


The above notation follows from the one level excitation model defined 


in Appendix D. 
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APPENDIX F 

Diffusion Approximation for Electron Leakage 


Instead of f(r,u,t) in equation (II^-l), consider an equivalent form 

■— A A —A A >■ . A 

f(r,ft,E,t) where f(r,ft,E,t) d dEdft is the number of electrons at r in d r 

A A 

whose flight direction ft is in the differential solid angle dft around ft 
and whose energy is between E and E + dE. With the differential electron 
flux defined as: 

<Kt, A.E.-t'USLde = ((t Sl.E/Ovf ci-TL dE (p.jj 

The Boltzmann equation [Equation (II-l)], in the absence of an external 
force field, becomes 


iA = -di-NT $ Cv - Z t CE\ §C*,il,e, t) 

-4- [[ n'e’-tUe'dfl.' + sCf.il ,e ,t) 

JJ ' (F-2) 


where we have replaced the collision term by the last three terms. In 
equation (F-2) Z^. (E) is the total collision cross-section at energy E and 

A A 

Z (E’ -*■ E, ft' - ft) is the cross-section for a process in which an electron 
with energy E' and the direction ft' is scattered into the element of solid 

A 

angle dft around ft and the energy interval dE at E. The term S(r,ft,E,t) 
accounts for other primary or secondary electron sources. We assume rota- 
tional symmetry around a distribution axis arbitrarily defined as the x-axis. 
Then the steady state form of equation (F-2) becomes: 




- 3 ) 
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We expand the angular dependence in Legendre polynomials, i.e.; 




(F-4) 


where 


SO.il, =Z iZiU s„(x.e;R,(^s'> 


/I /i 

Ge>s ©o — £2- m fL 


(F-5) 


(F-6 ) 


By the addition theorem for Legendre polynomials, Pn(Le>s©,P) can be 
expressed in terms of pJ^Ccose^ and Cgos©'} where the are the 

associated Legendre Polynomials. Performing the angular integral of 
equation (F-3) and using conventional orthogonal relations one arrives at 

t * ie approximation : 

fn + n = 

XM+W 0+1 'ax 

Czh+i) f 22 c x. e'-» g)$ Cx £' ) ^6.' -+ C^ n+ 0 Cx,e) 


n ~ °^M. 


(F-7) 


Truncating the equations with n = 1, we obtain the P^ approximation or the 


diffusion equation. For a homogeneous medium, this becomes* 


[.62 ] . 


-Me^^Cpe) + E-tCe) - [ r 4 U’^E)?( f.e'jdE 1 -t- 


(F-8) 


where 


cfOr. e) - | )«UX- 


( r , ^ s , 



■ u 
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X^Ce'^e)- j r s (E , -*e y Q. 4 ')clSl u 

4n 

and 

T>(e} = 1 — 

S(r t (E)'I s Cne.) (F-9) 

In the calculation, we have made the simplying assumption of isotropic 

scattering so the averaged cosine of the scattering angle vanishes. This 

form for the diffusion coefficient is only valid when the electron motion 

is not influenced by electric field forces. This restricts it to large 

electron kinetic energies and a small degree ionization. Furthermore, the 

diffusion equation (F-5) carries the usual restrictions for the validity of 

Fick's Law 1 J , one being that the size of the system be large relative to 


the mean free path of the electrons. For a 1-cm radius tube, this restricts 
the gas pressure to be, roughly, above one torr. 
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APPENDIX G 

Energy Deposition Rate and Primary 5 - Ray Energy Distribution 
by a Two Region Model in the Plasmas Created by Nuclear Radiation 


Using a two-region model, J. Guyot^^and q. Miley^^have derived 

expression for the rate of energy deposition by heavy-charged particles in 

gaseous media. In the experiments performed by J. Guyot^^and T. Ganley^^, 

laser tubes coated with B-10 were placed in the University of Illinois ■ TRIGA 

Reactor. The pulsed neutron flux falls upon the boron coating inducing the 

nuclear reaction n«f- B — * Li-t-«. . The heavy-charged particles traverse the 

and deposite their energy as they slow down. 

The geometry of the system and the heavy-charged particle reactions are 

shown in figure (G-l) where a slab geometry is used to approximate the cylinder 

til 

For 0^x^x cut , the rate of energy deposition R^(x.T q ) at x for the i type of 
heavy-charged particle with initial energy T q is: 


n+r tr>- ^°.c >] 




(G-l) 


and for x ^ ^ x s A (T ) , 
cut ' g o 


ft ; t ^) 1 

^ *- AcCla') -1 


where in (G-l), (G-2) ; S(T q ) = source rate of heavy-charged particle 

x = distance from the surface of coating 

x cut = E A(r.)) 


(G-2) 
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B- 10 GAS 



<L 


Figure G-l. Slab Geometry of the Source and Gajeous Medium 
Equivalent to the Actual Cylindrical Geometry 
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A (T q ) = range of heavy charged particle at T q in the Boron Coating 

= range of heavy charged particle at T q in the gas 

r = thickness of Boron region 

J a U- ^ J 

h = Cl- 

= ( VaU) - VA^CnO 

n = a fitting parameter for the slowing down and energy degradation. 

\ "l ^ 

For n^o , ■=• Li^Cl- 1 ^^ - 17T\ J ^ 

The following parameters are used in the calculation: 

T q (MeV) Ag(T o )cm A(T o )cm n 

a 1.459 2738.7/p 3.7 x 10 " 4 0 

Li 0.855 1240/p lO -4 -0.25 

3 -2-1 

and S(T q ) = 450 R 0/cm sec, where 0 is the neutron flux (cm -sec ), 

R is the fraction of the type of heavy charged particles. For the production 

of 1.459 MeV a and 0.855 MeV Li, R = 92%. The coating thickness is 0.4 mg/ cm". 

Figure (G-2) shows the rate of energy deposition versus gas pressure at 

the center of the tube. The tube diameter is 2.54 cm and the neutron flux 

13 -2 -1 

used is 0 = 5 x 10 cm sec . The rise in energy deposition rate as 
pressure increases is because of the increased slowing down efficiency as 
gas pressure is increased. This reaches a maximum at 'v 550 torr because 
fewer heavy-charged particles can reach at the centerline as pressure goes 
higher. 

The heavy-charged particle flux spectrum F^(x,T,T^) for the i*-^ type 
heavy charged particle with initial kinetic energy T q at x with kinetic 
energy T is given by: (133 ff. ref [64]) 



(KeV/cm 3 sec) per Neutron F 
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where 


e = T/T q , z = x/Xg(T o ), S = t/X c (T q ) 


e = (1-z) 
max v 


1/n+l 


£ . = (1-S-Z) 

cnt 1 


z . = (1-S) 

cnt 


1/n+l 


The primary 6-ray electrons are created by the passage of this flux 

through the gas. The spatial and energy distribution of these primary 

■f v» 

6-electrons is related to the i La type heavy charged particle flux by: 

T- 

" (x.T.WZ.CT, Wl 

(G-4) 


n.CV,H-0)=5 F; Cx,T,To)Z. t CT, E)<i7 

where T t ^ res ^ is the minimum kinetic energy for the i^h type heavy- 
charged particle to ionize and E^(E,T) is the probability per unit flux 
that a particle with energy between T and T + dT will lose an amount of 
kinetic energy between E and E + dE through ionization. Finally, U is the 
ionization potential of the gas atoms. Employing the formulation derived 
by M. Gryzinski for £^(T,E) and T t h res ^ (equations (6) and (8) of reference 
[23] ) , n^(x,E) has been solved for x = centerline of the tube at 1.27 cm and 
for 1.495 MeV a and 0.855 MeV Li ions which account for 92% of the heavy 
charged particles produced in the neutron irradiation of boron. Figure (G-3' 
shows the primary 6-electron production rate, per unit neutron flux, as a 
function of energy in 10 torr helium. 


1 

mi*.* 
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The total electron density is determined by a balance of the source 
due to the heavy-charged particles and the losses due to recombination and 
leakage. The main body of electrons is assumed to be thermalized and the 
losses are mainly due to the recombination and diffusion leakage of thermal 
electrons. With this assumption, the thermal diffusion coefficient is taken 
to be the ambipolar diffusion coefficient. For a 1/U recombination cross- 
section, with only one species of ions present, the steady state electron 



where D & = ambipolar diffusion coefficient 

R = rate of energy deposition by heavy charged particle slowing down 
- ™-'‘a R ~="recombination coefficient 

w = energy deposition for each ion-pair production 
A = first fundamental mode characteristic length of system 
However, if molecular ion formation by three body collision is included, 
electron recombination is predominantly with the molecular ion. Then the 
electron density becomes: 


- 2 . 0 . 


— : - -v- -l -"i ' 

A*- 1 V itfec*. 


(G-6) 


s + n 

A 

-3 -1 

S g = R/w = rate of electron production/volume [cm -sec ] 


D^, = ambipolar diffusion coefficients for the atomic ions and 

2 

molecular ions, respectively [cm /sec] 

rj = coefficient of molecular ion formation through three body 
collisions [1/sec] 



Iron Dsnslty (crrf^xlC 12 



Pressure (torr) 

Figure G-4 Electron Density vs Pressure at the-Center of Tube for Tube Radii = 0.635, 1.27 
cm and a Neutron Flux of 5 x 10 J cm" -sec 
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Pressure (torr) 

Figure G-5. Electron Density vs. Pressure a^the Center of Tube for Tube Radii = 0.635, 1.27 
cm and a Neutron Flux of 5 x 1Q 1 ^ cm" -sec" 
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L. 


Equation (G-6) has been evaluated using the coefficients given in references 
[ 49 ]* [ 50 ] and [ 63 ] and the resulting electron density vs. pressure for 
several tube sizes and neutron flux densities is illustrated in Figures 
(G~4 ) - (G-6 ) • The densities shown are for the centerline of the tube. The 
cross-over of the curves at about 30 torr is due to the balance between 
leakage and energy deposition rate. The peak electron density occurs at 
about 550 torr because higher pressures reduce the number of heavy-charged 
particles reaching the centerline. 

The addition of a distributed primary source S(E) in the low energy 
region makes equation (III-10) inhomogeneous; namely: 


(G-7) 

A series solution $(y) for equation (III-10) was given earlier in equation 
(III-15). [A second independent homogeneous solution , (y) , given by 

De Sabre (equation (8) of Ref. 65) was discarded by Horwitz et all'll because 
it gives negative slowing down density at y = 0.] A particular solution to 
equation (G-7) can be found by the method of variation of parameters from 
the two homogeneous solutions. This method gives: 


where 




(G-8) 


Equation (G-8) has been evaluated numerically, the particular solution 
added to the homogeneous solution, 4^(y), and the entire distribution nor- 
malized. Notice that as y -*■ 0, ^(y) = 0 and ^(y) again has a zero 
coefficient. Results are presented in figures fc- 7 ) and (b_g ) along with 
a plot of the distribution for the homogeneous case. It is seen that for 


Normalized Distribution 




{JTn t)/n 





lew energies where the electron population is high, the deviation is 
insignificant. For higher energies, where the population is already 
small,- the difference is about a factor of three. 
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APPENDIX H 

Series Solution of the Thermalization Equation With 
Recombination and Leakage 


-With the notation of Chapter III, the thermalization equation, 
the second moment of AE, becomes: 

-Dte'W'ipie.) + ="5 XL S ( e) [( z kT - B ^ -v 


with = ^Ce)MCE^ 


d^Ct) _ $ Ce) i + _L A_ 

de KCE) L hi E i MCE) dE 


d tc K e) 


i 

E 2 - 


tCe} + 

HLE) 




i_ afw 

HCE) dE 



$(E)1 J— 4_ §CE) 

Fci) J + de 2 - 


From equations (H-l), (H-2) and (H-3) , we obtain 


'XeO^a(E) ^Xvie)iCE) =7 z s Ce)L§Ce) + e^6) + fAT^# C e) . 


To change this result into the variable y, where y = /^Tb T , 

note that : 


we 


d. _ J L d 

de ~ 2.KJ dj 3>i 



and with these relations, equation (H-4) may be written as: 


y0' a ' t Ci-y-3v-A]$(p=-o 


including 


(H— 1 ) 


(H-2) 


(H-3) 


(H-4) 


(H-5) 


(H-6) 


where g and A are defined in the text. The solution for the special case of 
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no leakage or recombination, with the boundary condition $(y) = 0 at y = 0 
2 -y 2 

is c y e ; i.e. a Maxwellian distribution. We can use this to find a 

correction factor S(y) applicable to the general case. If one lets 4>(y) = 
2 -y 2 

ye 7 S(y), the differentials of <J> (y) with respect to y are: 




(H-7) 


$“0j) - Cso/H*- i o H a + 4^ j + s'cy)(^y- ( H -s) 

4 yM + s"<?) y 2 - 3 

The differential equation for the series S(y) from equations (H-6) to (H-8) 


is : 


V 5 ( y) ■+ 1 2 sUy) - sc/)[4 = o 

«=»o 

Writing out equation (H-9) in terms of the series S(y) = ^ZA 0 V / 

»nzo • 

we obtain: 

The coefficients of (H-10) are seen to be: 

A. = — [*■(*-!■->*/•) A,. 


(H-9) 


(H-10) 


(H-ll) 


With the normalization such that /” $(y) dy = 1 , the first five coefficients 
of are: 

A. -Vji , A, = § a. , + 

A s * a? r* r <- 1 ) */-> -+ w* -+ , 

and ] . < H - I2 > 

? ? 

Collecting terms of order 1, A, g, A“, g“ and Ag; for A, g << 1; we 


can write S(y) as: 
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SOp = » -V + f 1 ^) * 4 ^) + ^X(yj](H- 13 ) 

The series with first order in g appears only with even powers of y; 

*!» 0| = 51 b \> - b ‘MjnTo 1 = ^ 

' n=i zr > t m z(n-0 

-pf - if v'— • 


i .e . 


(H-13) 


The series with second order in g appears only after y and appears only in 
even powers of y: 

V'P = ^ ^*Cn+i) ^ (H-14) 

where j _ — ! P^z/PcL -t [) 2n 1 

d Hn+t) ~ (tn+2 )(za-r<D L J *” J 

and d 4 =l/4-6-8 

This gives 4 t = >4-«r^/4..6-TJ ^ I/6-&1 , d* - A** £!A-S? CZl + Ag)] 

The series with first and second order in A have been developed by 
Cohen and they are given as: 


•X u z ^i 

Aty)- Z - 

' »--• ( 2 o-m) p (n + r A) 


(H-15) 


and 


< 3 “ , T ->1 


m««)i [f. 


(H-16) 


The series with terms in order of Ag starts with y and appears only for 
odd powers of y. The coefficients of this series are: 

p» + 1 


=m-V)2: e„„H 


° - 1 


where 


p _ p | 2. C 2r\-0 

~ ^2-vi 


L( zn+ 


. e < = i - 


(H- 17) 
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The first few terms of c 2^ t( are; ^3 - i: + "g 

^s- = ic4- + ■?) C4r s )(rjij) , e 7 - i (-f + i ^ '^7 V - ) 

These series converge rather rapidly for y < 5. For larger values of 

y, the series is given as an inverse power polynomial of y. As it converges 

rapidly for large y, the series need not be broken into orders of g and A. 

r6 3 1 

De Sabrino et al give two independent asymptotic series for large y: 

top = r 1/i £c„ y" 

C. - “V*. , c. OAy 

C _ 4C,., (<\k t n-t-X +"-4) 

11 in 2n 


with 


(H-18) 


and 


with 




cp 


y-V4 & P 


Eaj 

/ 


-n 
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(H-19) 


The series <Hy) is damped by the Maxwellian factor; <5(y)— y Z <61 ^ (y) 

and it becomes negligibly small because of the exponential factor. The 
asymptotic solution for equation (H-6) is, therefore: 


$>0j) - 


= NJ 


-If 


v ‘£ 

n-o 


2 _a„y 


-71 


(H-20) 


Effect of Electric Field on High Energy Electrons 
During Slowing Down 


Consider an electron above the ionization energy of helium in an electric 

f <- ; 1 

field E eV/cm. The drift velocity u , due to the field is (pg. 122 of ref. 


hl e =e E 

c dt 




u-n 


where x is the mean collision time. The electrons drift distance in the 
direction of the field in this time is ^ ^ , and they acquire an energy 

AT given by: 


rvi . 


AT = SE 7 ! 1 _ 


e E 


21 7 - 


Cl-2) 


where E^ is the macroscopic inelastic cross -section and c the random speed 

of the electron. Consider two extreme cases of 1 Kev and 30 eV electrons in 

helium at 1 Torr with an electric field of 10 V-cm 1 . For 1 Kev, E - 0.5 cm 

(pg. 63 of ref. ^ ) and AT = 0.2 eV. For a 30 eV electron, E = 3 cm * and 

AT = 0.3 eV. The loss of energy due to an inelastic collision, however, is 

larger than 20 eV. For higher pressures, the gain in energy due to the drift 

2 2 

in electric field is even less because AT ~ 1/E ' 1/p . After an inelastic 

collision, the electron emerges in random directions. The drift in the elec- 
tric field has to be started anew. As illustrated above, the gain in energy 
from the field between two inelastic collisions, for both of those extreme 
energies, are less than 1 eV. It is concluded that the energy gain in the 
field can be neglected compared to the energy losses for electrons with 
energies above the ionization energy. This, however, is limited to moderate 
electric fields of 10 V-cm * or less. 
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APPENDIX J 

Discussion of the Equation for the Distrubtion Function 

in an Electric Field 

Various authors ^ 1 J ’ f 2 ] > [ 3 ] use d different forms of the equations 

for the two term expansion f 0 (u) and f ^ (u) of the distribution function in 


an electric field. Chapman and Cowling 
f Q and f^ in velocity space as: 


[67] 


give the coupled equation for 


~ 'S . ALA, ,‘rrr 


(J-D 


(J-2) 


with 


f 0 ^ = V 0 Cc) 

£,C C > - 

c. 

0<O ^ C /-^Cc) 

where ^Cc*), ^ and "$_Cc) are defined in chapter IV and the 

leakage term has been neglected. Without recombination Equation (J-l) 


becomes : 




(J-3) 


and the LHS of equation (.J-2) is: 

- AilcH4vbc?)] 


(J-4) 


With v (c) expressed as E N.cQ. (c) where j denotes the type of molecules, 

C i J J 

J r 31 

N. the density and the momentum transfer cross section; Nigham 1 J arrives 

at a set of coupled equations: 
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. A 

"iroo 1 (E 


» ' Xu ( °0 ~ [t^) c 


" 4j ^ ( ^ * "21 M, i CoJof^Oj (J-6) 

J 

Combining (J-S) and (J-6) and changing the variable u to UJ - , with 

d _ d- > Nigham obtains the following working equation: 

Ctu e 

”3 aw /nty (? ej Ml = (j-7) 

However, if recombination is included, equation (J-l) becomes 

r ico “ilS 

' ,C * " L eycKc) i ■*<- CJ-8) 

r h 71 

and equation (J-2) is given by Chapman 1 as: 


3 V m e J 1 I'AKc-) •'A i.Cc) J / 


(J-9) 


where unlike in Nigham' s case, ionization sources, recombination and the kT 
term have now been included. Combining equations (J-8) and (J-9) and 
assuming 8 £■ (c) << c, we find: 

r <- 

*- hr*, jko l& Cc) - £,CO = \ U-/2?JcMc fJ _ 9} 

MjlW J *0 MKO yj 

Then, after a formal integration, equation (IV-6) of chapter IV is obtained. 


121 


LIST OF REFERENCES 


1. T. Holstein, ’’Energy Distribution of Electrons in High Frequency Gas 
Discharge,” Phys. Rev. 70 , 367 (1949). 

2. M. J. Druyvestyen and F. M. Penning, "The Mechanism of Electrical 
Discharges in Gases of Low Pressure," Rev. Mod. Phys., 12 _, 87, (1940). 

3. W. L. Nighan, "Electron Energy Distributions and Collision Rates in 
Electrically Excited N 2 , CO and CO 2 ," Phys. Rev., 2A , 1989, (1970). 

4. V. Fano, "Degradation and Range Staggling of High-En rgy Radiation," 
Phys. Rev. 92_, 328,(1953). 

5. R. T. McGinnies, "Energy Spectrum Resulting From Electron Slowing 
Down," Cir. 597, National Bureau of Standards, Washington, D.C., (1959). 


6. G. R. Russell, "Feasibility of a Nuclear Laser Excited by Fission 
Fragments Produced in a Pulsed Nuclear Reactor," Trans. Symposium on 
Res. on Uranium Plasmas and Their Tech. Appl., University of Florida, 
Gainesville, (1970). 

7. D. R. Rees, C. B. Leffert and D. J. Rose, "Electron Density in Mixed 
Gas Plasmas Created by Fission Fragments," J. Appl. Phys., 4£, 1884, 
(1969). 

'8. K. Thom and R. Schneider, "Nuclear Pumped Lasers," AIAA, 10_, 400, (1972). 
Also J. C. Guyot , G. H. Miley, J. T. Verdeyen, and T. Ganley, "On Gas 
Pumping via Nuclear Radiations," in Research on Uranium Plasmas and 
Their Technological Applications, NASA SP-236, NASA, Washington, D.C., 
(1970) . 

9. N. Krascella, "Analytical Study of the Spectral Radiant Flux Emitted 
From the Fuel Region of a Nuclear Light Bulb Engine," U.A.R.L. Report 
J-910904-1 , (1970). 

10. E. P. Muntz, "The Electron Beam Fluorescence Technique," NATO Report, 
AD-692-860, (1968). 

11. J. T. Ganley, J. Verdeyen and G. Miley, "Enhancement of Co ? Laser Power 
and Efficiency by Neutron Irradiation," App. Phys. Lett., 18 , 568, (1971). 

12. J. C. Guyot, G. Miley and J. Verdeyen, "Metastable Densities in Noble- 
Gas Plasmas Created bv Nuclear Radiations," J. App. Phys., 42 , 5379, 
(1971). 

13. J. C. Guyot, "Measurements of Atomic Metastable Densities in Noble Gas 
Plasmas Created by Nuclear Irradiations," Ph.D. Thesis, University of 
Illinois, (1970). 

14. B. Wang, G. Miley, "Monte Carlo Simulation of Radiation-Induced Plasmas" 

Trans. Am. Nucl. Soc. 15, No. 2, (1972). ’ 


122 


15. B. Wang, "Monte Carlo Simulation of Non-linear Radiation Induced 
Plasma," Ph.D. Thesis, Dept, of Computer Science, Univ. of Illinois, 
Urbana, Illinois, (1972). 

16. L. V. Spencer and V. Fano, "Energy Spectrum Resulting From Electron 
Slowing Down," Phys . Rev. 9_3, 1172 (1954). 

17. L. F. Miller, R. E. Faw, "Effects of Radiation Quality on the Radiolysis 
of Water," Sp. Report #66, Kansas State Univ., Manhattan, Kansas, (1966). 

18. C. Moller,"Uber den Stoss zweier Teilchen unter BerUcksichtigung der 
Retardation der Krafte," Zeit F. Phys. 70_, 786, (1931). 

19. L. Vriens , "Binary-Encounter Electron-Atom Collision Theory," Phys. 

Rev. 141, 88, (1966) . 


20. L. Vriens, "Electron Exchange in Binary Encounter Collision Theory," 
Proc. Phys. Soc. 89_, 13, (1966). 

21. Von F. Lenz, "Zur Streuung mittelschneller Elektronen in kleinste 
Winkel," Z. Naturforschung 9A , 185, (1954). 

22. M. Gryzinski, "Classical Theory of Electronic and Ionic Inelastic 
Collision," Phys. Rev. 115, 374,(1959). 

23. M. Gryzinski, "Two Particle Collisions. I-General Relations for 
Collisions in the Laboratory System. II-Coulomb Collisions in the 
Laboratory System of Coordinates," Phys. Rev. 138 , A305-335, (1965). 

Also, "Classical Theory of Atomic Collisions. I-Theory of Inelastic 

Collisions," Phys. Rev. 138 , A336-358, (1965). 


24. E. Gerjouy, "Cross-Section for Energy Transfer Between Two Moving 
Particles," Phys. Rev. 148, 54,(1966). 


25. R. C. Stabler, "Classical Impulse Approximation for Inelastic Electron- 
Atom Collisions," Phys. Rev. 133, A1268, (1964). 


26. V. I. Ochkur and A. M. Petrunkim, "The Classical Calculation of the 
Probabilities of Excitation and Ionization of Atom by Electron Impact," 
Opt. Spectry. (USSR) T4, 245, (1963). 

27. A. Burgess, Proceeding of the Third International Conference on Elec- 
tronic and Atomic Collisions, North-Holland Publishing Co., Amsterdam, 
p. 237, (1964). 

28. C. B. Opal, W. K. Peterson and E. C. Beaty, "Measurements of Secondary- 
Electron Spectra Produced by Electron Impact Ionization of a Number 

of Simple Gases," J. of Chem. Phys. 55 , 4100, (1971). 

29. D. Rapp, P. Englander-Golden, "Total Cross Sections for Ionization and 
Attachment in Gases by Electron Impact. I. Positive Ionization," 

J. Chem. Phys., 4 3 , 1464, (J965). 

30. Y. Kim and M. Inokuti, "Total Cross-Section for Inelastic Scattering of 
Charged Particles by Atoms and Molecules. V Evaluation to the Next Order 
Beyond the Lethe Asymtote , " Phys. Rev. 3, 3, (1971). 


123 


31. Y. Kim and M. Inokuti, "Generalized Oscillator Strengths of the Helium 
Atom.. I," Phys. Rev. 175 , 176,(1968). 


32. L. Vriens, J. A. Simpson and S. R„ Mielszarck, Test of Born Approxi- 
mation: Differential and Total 2 J S, 2^P and 2^S Cross Sections for 

Excitation of He by 100 to 400 eV Electrons," Phys. Rev. 165 , 7, (1968). 


33. R. St. John, F. L. Miller, and C. C. Lin, "Absolute Electron Excitation 
Cross-Sections of He," Phys. Rev. .134, 888, (1964). 


34. H. S. W. Massey and B. L. Moiseiwitsch, "The Application of Variations 
Methods to Atomic Scattering Problems, IV. The Excitation of the 2^S 
and 2^S States of Helium by Electron Impact," Proc. Roy. Soc. A227, 

38, (1955). 

35. H. Maier-Leibnitz , " Aus b eut erne ssuh gen beim Stoss langamer Elektronen 
mit Edelgasatomen," Zeit. f. Physik, 95 , 499, (1935). 

36. T. Itoh and T. Musha, "Monte Carlo Calculations of Motion of Electrons 
in Helium," J. Phys. Soc. of Japan, 15 , 1675, (1960). 

37. W. Jesse and J. Sadanskis, "Ionization in Pure Gas and the Average 
Energy to Make an Ion Pair for Alpha and Beta Particles," Phys. Rev. 97., 
1665, (1955). 

38. L. R. Peterson and A.E.S. Green, "The Relation Between Ionization 
Yields, Cross Sections and Loss Functions," Proc. Phys. Soc. 1_, 1131, 
(1968). 

39. N. Mott and H.S.W. Massey, The Theory of Atomic Collisions , 3rd ed., 
Oxford U. Press, London, pg. 455, (1965). 

40. R. Labahm, and J. Callaway, "Distortion Effects in the Elastic Scattering 
of 100-400-eV Electrons From Helium," Phys. Rev. 180 , 91, (1969). 

Also L. Vriens, C. Kuyatt and S. Mielczarek, "Test of Born Approximation; 
Differential and Total Cross Sections for Elastic Scattering of 100- to 
400-eV Electrons by Helium," Phys. Rev. 170, 163, (1968). 


41. D. Golden and H. Bandel, "Absolute Total Electron-Helium Scattering 

Cross-Sections for Low Electron Energies," Phys. Rev. 158, 14, (1965). 


42. C. Normand, "The Absorption Coefficient for Slow Electrons in Gases," 
Phys. Rev. 35, 1217, (1936). 

43. H. Hurwitz, M. S. Nelkin and G. J. Habetler, "Neutron Thermalization , 
I. Heavy Gas Modulator," Nucl. Sci. Eng. 1_, 280, (1956). 

44. K. Beckurts and K. Wirtz , Neutron Physics, Springer=Verlag Press, 

New York, pg. SS, (1964). 

45. T. F. O'Malley, "Extrapolation of Electron-Rare Gas Atom Cross Section 
to Zero Energy," Phys. Rev.' 150, 1020, (1968). 


L. 


124 


46. D. Rose and M. Clark, Jr., Plasmas and Controlled Fusion , The MIT Press, 
Cambridge, Mass., pg 165ff, (1961). 

47. L. Spitzer, Jr., Physics of Fully Ionized Gases , Interscience Publ., 

New York, 2nd ed., pg 127ff, (1964). 

48. E. H. Holt and R. E. Haskell, Plasma Dynamic s, The Macmillan Co., 

New York, pg 269ff, (1965). 

49. H. J. Oskam and V. R. Mittelstadt, "Recombination Coefficient of 
Molecular Rare-Gas Ions," Phys. Rev. 132 , 1445, (1963). 


50. A. V. Phelps and S. C. Brown, "Positive Ions in the Aterglow of a Low 
Pressure Helium Discharge," Phys. Rev. 86 , 102, (1952). 

51. A. von Engel, Ionized Gases , Oxford Press, 2nd ed., pg 164, (1965). 

52. A. E. Heyen and T. J. Lewis, "Electron Energy Distribution Functions 

and Transport Coefficients for Rare Gases," Proc. Roy. Soc. 271, A531, 
(1963) . 

53. W. P. Allis, Handbuch der Physik , Vol. XXI, 304, Springer, Berlin, (1956). 

54. von J. A. Smit, "Berechnung der Geschwindigkeits-Verteilung der Elek- 
tronen bei Gasentladungen in Helium," Physica 3_, 543, (1936). 

55. Domenick Barbiere, "Energy Distribution, Drift Velocity and Temperature 
of Slow Electrons in Helium and Argon," Phys. Rev. 84_, 653, (1950). 

56. M. H. Hughes, "Electron Energy Distribution and Transport Coefficients 

in Helium and Neon," J. Phys. B: Atom. Molec. Phys. 3, 1544, (1970). 

57. J. T. Ganley, "Effect of a Heavy Particle External Ionization Source 

on Carbon Dioxide Laser Discharges," Ph.D. Thesis, Dept, of Electrical 
Engr., Univ. of Illinois, Urbana, 111., (1972). 

58. H. J. Oskam, "Microwave Investigation of Disintegrating Gaseous Discharge 
Plasmas," Philips Res. Report, 13, 432, (1958). 

59. A. K. Bhattacharya , Microwave Investigation of Plasmas Produced in a 
Nuclear Reactor, Ph.D. Thesis, Nuclear Engineering Program, Univ. of 
Illinois, Urbana, 111., (1966). 

60. N. Mott, "The Collision Between Two Electrons, " Proc. Roy. Soc. A-126, 

259 (1930). 

61. E. A. Hylleraas, "Neue Berechnung der Energie des Heliums im Grundzu stande , 
sowie des tiefsten Terms von Ortho-Helium," Zeit f. Phy. 54_, 347, (.1929). 

G. Bell and S. Glasstone, Nuclear Reactor Theory , Van Nostrand Reinhold 
Press, New York, pg 144, (1970). 


62 . 



12S 


63. J. C. Guyot, G. Miiey and J. T. Verdeyen, "Application of a Two-Region 
Heavy-Charged Particle Transport Model to Noble-Gas Plasmas Induced by 
Nuclear Radiations," Nucl. Sci. § Engr. 48 , 373, (1972). 

64. G. H. Miiey, Direct Energy Conversion of Nuclear Energy, Am. Nuc. Soc., 
Hinsdale, 111., pg 113ff, (1971). 

65. Luis De Sobrino and M. Clark, Jr., "A Study of Wilkins Equations," 

Nuc. Sci. Eng., 10, 388, (1961). 

66. E. R. Cohen, "The Neutron Velocity Spectrum in a Heavy Modulator," 

Nuc. Sci. Eng., 2_, 227, (1957). 

67. S. Chapman and T. Cowling, The Mathmetical Theory of Non-Uniform Gases , 
2nd. ed., Cambridge University Press, London, pp.356, (1960). 



126 


Ronnie Hung-Kie Lo was born in 

He received the B.S. degree in Applied Mathematics and Engineering 
Physics with honor from the University of Wisconsin in 1966. He received 
the M.S. degree in Nuclear Engineering from the University of Wisconsin 
in 1968. He joined the Picker Corporation and served as a R. § D. 
engineer for one and a half years. 

In September 1969, he came to the University of Illinois to study 
for a Ph.D. degree in the Nuclear Engineering Program. 

As an undergraduate, Mr. Lo was awarded the Non-resident Scholarship 
from 1963-66. In graduate school he was awarded the University Fellowship 
and Research Assistantships . Mr. Lo has co-authored with Professor Miley 
an article in the Transaction of the American Nuclear Society , 18th Annual 
Meeting, June, 1972, and a paper accepted for presentation in the 25th 
Annual Gaseous Electronics Conference, University of Western Ontario, 
London, Ont., October, 1972. 






